Table 1. Census of CFHT HRCam Images 



Field 


Filter 


Exposures 


Airmass 


FWHM 


Central Field 


V 


4 x 600 sec 


1.19 - 1.35 


0'.'72 


Central Field 


R 


26 x 500 sec 


1.08 - 1.31 


0'.'73 


Central Field 


I 


4 x 500 sec 


1.35 - 1.48 


0'.'73 


East (Background) 


R 


12 x 600 sec 


1.01 - 1.35 


0'/58 



M87 Central Field (R < 23.2) 
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Table 2 . Photometry of Brightest Clusters 



X" Y" R V-R R-I X" Y" R V-R R-I 





-28.57 


-11.67 


19.03 


0.55 


0.74 


-36.81 


-3.01 


20.70 


0.47 


0.72 




3.70 


25.88 


19.30 


0.48 


0.70 


6.03 


15.56 


20.73 


0.49 


0.66 




10.59 


-36.66 


19.60 


0.60 


0.65 


5.28 


23.06 


20.73 


0.47 


0.75 




22.65 


-36.81 


19.70 


0.62 


0.66 


87.72 


-38.18 


20.74 


0.63 


0.60 


oo : 


-12.34 


-15.05 


19.83 


0.50 


0.63 


-24.04 


-42.40 


20.75 


0.48 


0.71 
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68.59 


-34.00 


19.84 


0.55 


0.65 


-58.04 


-19.33 


20.76 






21.54 


-6.39 


19.95 


0.46 


0.58 


-43.58 


-0.46 


20.77 


0.50 


0.77 


an 


12.84 


-34.67 


20.04 


0.54 


0.58 


-5.28 


-7.87 


20.77 






50.22 


-61.90 


20.09 


0.55 


0.47 


9.93 


-27.16 


20.79 


0.52 


0.86 


i— j , 


-56.79 


-75.42 


20.13 


0.71 


0.65 


35.23 


-9.19 


20.80 


0.41 


0.80 


<N ' 
Cn , 


14.15 


9.40 


20.18 


0.42 


0.78 


44.48 


-69.60 


20.80 


0.52 


0.47 


26.25 


-102.68 


20.26 


0.50 


0.46 


25.69 


11.83 


20.81 


0.46 


0.71 


7— 1 ■ 


21.28 


-82.58 


20.27 


0.45 


0.47 


-31.55 


-18.29 


20.81 


0.54 


0.79 


> ' 
1 1 

T— 1 . 


-0.69 


-103.44 


20.29 


0.49 


0.56 


-2.73 


-38.39 


20.81 




0.63 


-34.35 


-47.79 


20.29 


0.47 


0.60 


-4.62 


-77.45 


20.83 


0.53 


0.66 


(N ■ 


-10.82 


-60.20 


20.33 


0.51 


0.55 


-27.84 


-9.68 


20.83 


0.62 


0.61 


<s_»> , 
00 . 


-1.47 


-13.44 


20.35 


0.36 


0.81 


74.10 


-26.36 


20.84 


0.57 


0.75 


-45.77 


-4.28 


20.36 


0.61 


0.73 


-20.46 


-72.29 


20.85 


0.46 


0.62 


On ■ 


19.17 


-103.42 


20.37 


0.49 


0.52 


-64.80 


-40.48 


20.85 


0.40 


0.65 


^ : 


-8.35 


-92.65 


20.38 


0.52 


0.49 


35.53 


-41.50 


20.86 


0.51 


0.62 




-21.02 


-12.21 


20.39 


0.49 


0.67 


-0.27 


-119.11 


20.88 


0.47 


0.64 


6 ■ 


-59.17 


-25.51 


20.39 


0.97 


1.30 


14.28 


-19.57 


20.91 


0.45 


0.77 
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C/2 . 


-18.68 


-46.16 


20.41 


0.51 


0.64 


8.97 


-36.83 


20.92 
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0.74 


-0.90 


-33.62 


20.44 


0.47 


0.51 


26.99 


15.26 


20.92 


0.58 


0.67 
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-23.97 


-45.10 


20.44 


0.56 


0.60 


-39.97 


-36.56 


20.92 


0.47 


0.59 




-4.00 


-45.45 


20.46 


0.58 


0.71 


4.28 


27.58 


20.93 


0.45 


0.70 




-2.40 


-55.19 


20.48 


0.59 


0.69 


-67.76 


-55.01 


20.94 


0.45 


0.53 




-44.39 


-27.15 


20.49 


0.55 


0.64 


12.64 


-105.18 


20.95 


0.39 


0.48 




35.53 


-37.00 


20.50 


0.53 


0.73 


-0.10 


12.27 


20.97 


1.09 


0.15 




24.77 


-19.50 


20.50 


0.50 


0.58 


-29.90 


-34.75 


20.97 


0.50 


0.51 




21.55 


3.43 


20.51 


0.51 


0.62 


-36.85 


-90.51 


20.99 
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0.74 




-39.90 


-29.61 


20.53 
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6.38 
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Table 3 . Definition of Radial Zones 



Zone 


r (px) 


r (arcsec) 


A (arcmin 2 ) 


a 


m 


Ring 1 


70 - 150 


9.17- 


- 19.65 


0.242 


1.812 


23.11 


Ring 2 


150 - 300 


19.65 


- 39.30 


0.822 


1.465 


23.62 


Ring 3 


300 - 450 


39.30 


- 58.95 


0.856 


1.590 


24.16 


Ring 4 


450 - 1080 


58.95 


- 141.5 


2.855 


1.460 


24.28 


Bkgd 








4.96 


1.975 


25.19 
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Table 4. The M87 GCLF as a Function of Radius 



R 


Ring 


1 




Ring 


2 


Rin 


g3 


Ring 


4 


Bkgd 


19.2 


0.0 


± 


1 


7 


2.0 ± 


1.7 


0.0 


± 


1.7 


0.0 ± 


2.0 


0.0 


± 


1.0 


19.6 


0.0 


± 


1 


7 


1.0 ± 


1.7 


1.0 


± 


1.7 


0.0 ± 


2.0 


0.0 


± 


1.0 


20.0 


2.0 


± 


1 


7 


2.0 ± 


1.7 


0.0 


± 


1.7 


2.0 ± 


2.0 


0.0 


± 


1.0 


20.4 


1.0 


± 


1 


7 


6.1 ± 


3.0 


11.1 


± 


4.0 


9.1 ± 


3.8 


0.0 


± 


1.0 


20.8 


2.9 


± 


2 


1 


11.7 ± 


4.2 


5.5 


± 


3.0 


10.4 ± 


4.4 


3.0 


± 


1.7 


21.2 


7.1 


± 


3 


2 


14.1 ± 


4.6 


8.9 


zb 


3.7 


13.6 ± 


4.7 


1.0 


± 


1.0 


21.6 


4.1 


± 


2 


5 


11.3 ± 


4.1 


11.2 


zb 


4.1 


22.4 zb 


5.8 


0.0 


± 


1.0 


22.0 


5.1 


± 


2 


8 


15.1 ± 


4.9 


24.0 


± 


6.0 


27.9 ± 


6.7 


3.0 


± 


1.7 


22.4 


8.9 


± 


3 


8 


33.8 ± 


7.3 


19.2 


± 


5.4 


31.8 ± 


7.1 


2.0 


± 


1.4 


22.8 


15.0 


± 


5 


5 


29.1 zb 


6.9 


18.5 


± 


5.3 


63.9 ± 


10.0 


2.0 


zb 


1.4 


23.2 


22.2 


± 


9 


9 


35.8 ± 


8.4 


25.0 


± 


6.4 


44.1 zb 


8.8 


6.1 


zb 


2.5 


23.6 










29.0 zb 10.6 


21.0 


± 


6.4 


43.9 ± 


9.5 


11.3 


± 


3.4 


24.0 














17.1 


± 


7.7 


32.8 ± 


10.0 


18.8 


zb 


4.4 


24.4 




















21.5 ± 


15.4 


48.9 


zb 


7.3 



Total 68.4±13.5 191.1±19.5 162.5±17.3 323.3±28.3 96.1±10.3 
Rel.Area 0.049 0.166 0.173 0.576 1.000 
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Table 5. Cluster Metallicity Versus Radius 



r" 


N 


(V-I) 


a(V - I) 


[Fe/R] V -i 


ct [Fe/H] 


20.5 


14 


1.139 ±0.031 


0.122 


-0.34 ±0.18 


0.66 


28.4 


17 


1.164 ±0.036 


0.146 


-0.20 ±0.19 


0.79 


38.1 


18 


1.145 ±0.027 


0.155 


-0.31 ±0.20 


0.85 


47.3 


22 


1.139 ±0.038 


0.168 


-0.34 ±0.19 


0.91 


60.0 


21 


1.148 ±0.015 


0.163 


-0.29 ±0.20 


0.90 


78.5 
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1.045 ± 0.022 
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0.81 


99.5 
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-0.62 ±0.21 
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Table 6. Globular Cluster Data for Selected BCGs 



Galaxy Cluster S n S$ a d (km s" 1 ) ICM kT (keV) 



N0524 


CfA13 


-22 


1 


6.4 ± 1.4 


N1275 


A426 (Perseus) 


-23 


1 


4.3 ± 1.4 


N1399 


S373 (Fornax) 


-21 


7 


11.1 ±2.1 


N1407 


Eridanus 


-22 


1 


4.0 ± 1.3 


N3311 


A1060 (Hydra) 


-22 


4 


15.0 ±6.0 


N3379 


CfA68 (Leo) 


-20 


7 


1.3 ±0.7 


N4073 


MKW4 


-23 


1 


4±2 : 


N4486 


Virgo 


-22 


4 


14.1 ± 1.6 


N5128 


Centaurus 


-22 





2.6 ±0.6 


U9799 


A2052 


-23 


4 


20.7 ±6.5 


U9958 


A2107 


-23 


4 


12.2 ±5.6 



4.9 ± 1.1 


206 ±110(1) 


6.7 


(6) 


3.3 ±1.1 


1026 ± 80(2) 


6.7 


(6) 


10.6 ± 1.9 


310 ±30 (2) 


1.1 


(6) 


3.1 ± 1.0 


264 ± 30 (2) 


0.9 


(7) 


11.5 ±4.6 


610 ± 50 (2) 


3.9 


(6) 


1.0 ±0.5 


120 ± 70 (1,3) 






3.1 ± 1.5 : 


525 ± 60 (2) 


1.7 


(6) 


9.1 ± 1.3 


632 ± 40 (2) 


3.0 


(6) 


2.0 ±0.5 


90 ± 32 (4) 






15.9 ±5.0 


750 ± 180 (5) 


3.1 


(6) 


9.4 ±4.3 


622 ± 70 (2) 


4.2 


(6) 



References for Table 6. 

(1) Geller & Huchra 1983; (2) Fadda et al. 1996; (3) de Vaucouleurs 1975; (4) Hesser et al. 
1984; (5) Girardi et al. 1993; (6) Harris et al. 1995; (7) Davis & White 1996 
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ABSTRACT 



We have used the High Resolution Camera at the Canada-France-Hawaii 
Telescope to obtain VRI photometry of the globular clusters in the innermost 
140" of the M87 halo. The results are used to discuss several issues concerning 
the formation and evolution of globular cluster systems in supergiant ellipticals 
like M87. Our principal results £1X6 clS follows: (1) From our deep R— band 
photometry of the cluster population, we find no significant change in the 
globular cluster luminosity function (GCLF) with galactocentric radius, for 
cluster masses M> 1O 5 M . This result places constraints on current theoretical 
predictions of the rate of cluster evolution by tidal shocking and evaporation, 
indicating that the main effects of dynamical evolution may be only on 
lower-mass clusters ( < 1O 5 M ) that are below the faint limit of most current 
observations. (2) Combining our (V — /) color indices with other data in the 
literature, we derive the metallicity gradient and mean metallicity of the clusters 
from r = 9" out to r ~ 500" (~ 0.7 to 35 kpc). Within the core radius r c = 1' of 
the globular cluster system, the metallicity distribution is uniform, but at larger 
radii the mean metallicity declines steadily as Z/Z & ~ r ~ - 9 . (3) The various 
options for explaining the existence of high specific frequency galaxies like M87 
are evaluated. We argue that alternatives involving secondary evolution (such 
as the formation of many globular clusters during mergers, or the existence of 
a population of intergalactic globular clusters centered on the same location as 
the central elliptical) are capable of modifying the characteristics of the GCS in 
distinctive ways, but are not likely to be the primary explanation for high— Sn 
galaxies. (4) We offer a new explanation for the large observed Sn range 
among brightest cluster ellipticals (BCGs). It is suggested that these central 
supergiant E galaxies formed in an exceptionally turbulent or high-density 
environment which favored a very rapid initial star formation burst. As a result, 
a much higher than average fraction of the protogalactic gas was driven out 
in a galactic wind not long after the first globular clusters were formed, thus 
artificially boosting the specific frequency of the remaining galaxy. From a total 
sample of 30 BCGs, we derive empirical scaling relations which relate to this 
hypothesis. Our analysis favors the view that BCGs began forming at redshifts 
z>5, distinctly earlier than most other galaxies. 



Subject headings: Stellar Systems: Globular clusters; Normal galaxies 
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1. INTRODUCTION 

Although much is now known about the characteristics of globular cluster systems 
(GCSs) in large galaxies (e.g., |Harris 1991| ), the properties of globular clusters in the inner 



regions of these galaxies are still poorly understood, for the simple reason that such objects 
are extremely difficult to observe against the strong background light of their host galaxy. 
The innermost few kiloparsecs are, however, precisely the regions where the strongest traces 
of dynamical evolution of the clusters should be seen, and are thus of considerable interest 
as tests of the many dynamical models now in the literature (see, e.g., |Okazaki fc Tosa 1995 



Capriotti fc Mawley 19961 ; jGnedin fc Ustriker 1997] ; |Ustriker fc Gnedin 1997] ; |Vesperini 199?) 



for comprehensive recent discussions). Do the inner clusters exhibit the expected history 
of disruption through their spatial distribution, luminosity distribution (GCLF), orbital 
distribution function, or local specific frequency? What fraction of the globular clusters 
that might have initially formed in these inner regions survived successfully to the present 
day? At present there are few galaxies for which relevant GCS data exist to constrain the 
models strongly. 

From the theoretical side, an additional major problem is to decide which of the 
systemic properties of the GCS that we see today (spatial distribution, GCLF, and so on) 
are due to the formation epoch, and how much to subsequent dynamical evolution. There 
are now plausible theoretical and observational grounds to suppose that the formation 
mechanism for globular clusters in protogalaxies was largely independent of location in 
the halo ([Harris fc Pudritz 1994j [hereafter HP94]; [McLaughlin fc Pudritz 1996t [Elmegreen 



fc Efremov 1997). But as soon as clusters form, they fall prey to a variety of erosion 



mechanisms, especially tidal shocking, early mass loss driven by high-mass stars, and (in 
disk galaxies) disk shocking. In the innermost < 2 kpc, dynamical friction may also play a 
role in removing the most massive clusters (e.g., Capriotti fc Hawley 1996; Capuzzo-Dolcetta 



fe Tesseri 1997| ). All these mechanisms depend sensitively on galactocentric distance R 



',!<■ 



and the tidal field. Many theoretical discussions (cf. the references cited above) have raised 
the possibility that the inner-halo structure of the GCS might be completely dominated 
by dynamical evolution, i.e. that the great majority of the clusters originally present there 
might now be gone. For M87 particularly, several recent papers discuss the possibility that 
the entire massive nucleus might have built up from dissolved clusters (see |Capuzzo-Dolcetta 



1993| ; [McLaughlin 1995| ; |Capuzzo-Dolcetta fc Vignola 1997|) . An analysis specifically tuned 



to giant ellipticals such as M87 QMurali fc Weinberg 1997| ) also suggests that both the total 



population and mass spectral index of the GCS should suffer quite significant changes with 
time because the rates of dynamical erosion depend on galactocentric distance. 



On observational grounds, it is now clear that GCSs in large galaxies display spatial 
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distributions which typically have large core radii, typically a few kpc in size (e.g., Lauci 



fc Kormendy 1986| ; |Harris et al. 1991] ; |Grillmair et al. 1994a| ; |McLaughlin 1995| ; [Forbes 
et al. 199"6]). It has frequently been suggested that this feature might have resulted from 



wholesale cluster erosion at small R gc , a view which has been argued by the dynamicists 
since the early paper of Tremaine et al. (1975). Some recent continue to favor such an 
interpretation (e.g., Capuzzo-Dolcetta 1993| and |Capuzzo-Dolcetta fc Tesseri 1997] ), while 
other discussions favor considerably more modest effects of destruction (e.g. [McLaughlin 
1995|) . Clearly, the possibility that many clusters might have been destroyed in these inner 



regions does not require that a large initial population of clusters must have existed in the 
core regions of these galaxies; the question is still open, and the direct evidence that we 
are left with today is highly circumstantial at best. An alternate approach (though one 
which has not yet been explored in the same detail; see HP94) is that a similar result could 
have arisen from the details of the cluster formation epoch. For example, the strong tidal 
field could have inhibited the initial buildup of globular-like protoclusters, or early mass 
loss could have evaporated them rapidly in the first few ~ 10 s y of their history. Direct 
observational evidence of cluster formation in the inner regions of the merger remnant NGC 
7252 recently presented by Miller et al. (1997) is intriguingly consistent with such a view, 
and deserves much closer evaluation as well as theoretical modelling. 

A somewhat more direct piece of evidence that dynamical evolution has dominated the 
construction of the GCLF that we now observe - in essence, something closer to a "smoking 
gun" signature of erosive processes - would be a progressive change in the shape of the 
GCLF with galactocentric distance, i.e. the relative numbers of clusters of different mass. If 
the initial mass distribution of the clusters is at least roughly independent of galactocentric 
distance (e.g., HP94; |Elmegreen fc Efremov 19971 [McLaughlin fc Pudritz 1996| ; piFazio| 
1986 ; |Ostriker fc Gnedin 1997| ), then such a change can only have arisen from dynamical 
evolution. On the theoretical side, this is the approach taken by several other authors (e.g., 



|Capriotti fc Hawley 1996| , |Vesperini 1997| , or |Ostriker fc Gnedin 1997| ). We now need better 
and more extensive data to compare with these simulations. 

Within the Milky Way, the population of globular clusters at all R gc is quite well 
defined over the full range of luminosities, but is too small a sample to provide strong 
statistical tests (e.g., Gnedin fc Ostriker 1997 ; Kavelaars fc Hanes 1997 ). M87 itself is 
perhaps the best available testbed for these ideas, because of its simple gE structure, 
relatively nearby location, and enormous cluster population: statistically useful numbers of 
globular clusters can be traced in to within ~ 1 kpc of its center, and the GCLF can be 
studied to faint levels. In this paper, we present new photometric data for the inner-halo 
GCS of M87 specifically designed to address these issues. 
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The paper is cast in two major parts. In Sees. 2-4, we present our new photometric 
data for M87, the luminosity function of the inner-halo GCS, and the color and metallicity 
distributions. Then in Sec. 5, we take up the long-standing problem of understanding 
how M87 and other central-supergiant ellipticals can generate such exceptionally large 
populations of globular clusters in all parts of their halos (the "specific frequency problem"). 
After assessing various suggestions in the literature, we conclude that the basic GCS 
characteristics were most likely built in at birth, and suggest that they may have been 
influenced strongly by the early galactic wind in such galaxies. 

2. THE DATA 

2.1. Observations and Calibration 

The observations for this study were obtained at the Canada-France-Hawaii Telescope 
during 1993 March 28/29 and 29/30. The High Resolution Camera (McClure et ai. 1989|) 
was used at prime focus, with the red-sensitive SAIC1 CCD detector: this CCD (now 
decommissioned) had an 18 micron pixel size, a 1020 x 1020 format after bias subtraction 
and trim, readnoise rms 5.5e~, and gain 1.9 e~ per adu. With a scale of 0'.'131 per pixel in 
HRCam, the total field size was 2'.23 on a side, corresponding to 9.7 x (D/15) kpc if M87 
is at a distance of D Mpc. This size was comfortably enough to enclose the entire central 
region of the M87 GCS, which has a core radius r c ~ V ( McLaughlin 1995f) . 

For our deep GCLF observations, we used the R— band, which was the standard filter 
nearest the peak sensitivity of the SAIC1 CCD. The GCLF turnover (peak frequency of 
cluster numbers per unit magnitude) is at V ~ 23.7 ( |Whitmore et al. 1995| ), or R ~ 23.2 
for (V — R) ~ 0.5 (see below). To reach well past this limit while avoiding saturation from 
the bright background light of M87, we used a long string of relatively short exposures 
(26 x 500 seconds) with individual frames shifted along a sequence of 7 different positions 
separated by several arcseconds. The individual exposures were then re-registered and 
median-combined to construct a final deep image. This image, which we will refer to as 
the "central field," is shown in Fig. [l| with the smooth galaxy light subtracted. The field 
center is actually almost 1' due south of the M87 nucleus, a choice which was governed by 
the location of the best available nearby guide star to drive the HRCam tip-tilt correction 
during the exposures. 

Although the population of detected starlike images in the M87 central field is totally 
dominated by globular clusters, an accurate correction for field contamination requires a 
control field. We selected a separate background field located 15'09 E and 0'45 N of the 
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central field and took another series of deep R exposures there. At this projected distance 
from the galaxy center, the number of M87 globular clusters is < 2% of the level in the 
central field (see the wide-field GCS profiles of McLaughlin et al. 1994). 

Finally, we took shorter series of the central field in V (4 x 600 sec) and I (4 x 500 
sec) to measure the colors of the brighter globular clusters in the field. These were used 
to estimate the size of the GCS metallicity gradient in the core region. A summary of the 
observational material is given in Table 1. 

Master flat fields in all three filters were constructed from a combination of dome 
flats and twilight exposures, and the control field was used to check that these produced 
globally flat calibrations to well within 1%. Absolute calibration was accomplished with 
VRI exposures of standard-star fields in the clusters M67, M92, NGC 4147, and NGC 7006 
( |5child 1983 ; [Christian et al. 1985| ; Pavis 1990|) . Both nights were photometric, and images 



of all four standard fields were taken on each night. The zeropoints of the photometric 
scales in our M87 fields are judged to be uncertain by no more than ±0.02 mag in each 
filter, from the observed scatter of the bright-star aperture photometry used to transfer the 
calibrations. 



Couture et al. (1990) obtained BVI CCD photometry of ~ 270 clusters over a roughly 
similar radial region around M87. Unfortunately, we cannot compare our data directly with 
theirs, since their field (just to the north of the M87 nucleus) has almost no overlap with 
ours. The much wider-field study of [McLaughlin et al. (1994)| in the V band can, however, 



be used for comparison since our outer zone overlaps with their inner one. For 44 objects in 
common brighter than V ~ 22.2, we find A\/(MHH-HHM) = -0.035 ± 0.014, confirming 
that there are no major zeropoint errors. 



2.2. Photometry 

To facilitate the photometry of the faint starlike objects in the central field, we first 
removed the background light from M87 using the ellipse-fitting codes in STSDAS. Fig. 1 
shows the field after subtraction of the smooth and nearly circular isophotes. As is apparent 
from the figure, the model also successfully removed the bright semistellar nucleus of the 
galaxy, leaving only the nuclear jet, the large population of globular clusters, and a few 
faint background objects. 



The photometry was done with the normal suite of DAOPHOT II and ALLSTAR codes 
( [Stetson 1992 ) with a star-finding threshold set at ~ 3.5 times the standard deviation of 
the average sky noise across the frame. Two iterations of the object-finding and ALLSTAR 
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measurement were sufficient to capture essentially all the starlike objects on the frames, 
since the field is not crowded in any absolute sense. Fig. |2| shows the locations of the 
brighter (R < 23.2) objects on the frame, along with the field orientation: the +x— axis is 
(32?95 ± 0?5) eastward from due North. A ±15? sector around the nuclear jet was excluded 
from any further analysis. A few clearly nonstellar objects (faint background galaxies) were 
removed objectively from the detection lists by the use of the image moment r_ 2 defined in 
Harris et al. (1991). Exactly the same rejection criteria were employed for the background 
field. For the magnitude range 19.0 < R < 24.6, 892 starlike objects in total were measured 
in the central field and just 91 in the control field. Table 2 contains a partial listing of the 
final data, for the brightest objects: here, (X, Y) are in arcseconds relative to the center of 
M87, with X increasing eastward and Y northward. An electronic file of the complete data 
may be obtained from WEH on request. 

Since the pixel-to-pixel sky noise depends dramatically on radius from the galaxy 
center, so does the completeness of detection at any given magnitude. (Object crowding 
is not important at any radius, so the completeness / is determined almost purely by the 
background noise). Extensive artificial-star tests were carried out to determine / and the 
internal precision of the photometry: scaled point-spread functions (psfs) were added to 
the frame, typically in groups of 200 — 300 over the entire magnitude range of interest, and 
then subjected to exactly the same measurement procedure as the original frame. If N(m) 
denotes the number of artificial stars inserted at magnitude m and N meas is the number 
successfully recovered (at any measured magnitude), then we define f(m) = N meas /N(m). 

After preliminary trials with various radial binnings, we divided the central field 
into four annuli as listed in Table 3 and marked in Fig. 2. The very innermost circle 
(r < 9" ~ 0.65 kpc [D/15]) was eliminated from further consideration because of its 
overwhelmingly high background noise, which prevented any useful photometry. The 
results for /, in steps of 0.5 mag, are shown in Fig. |3|. Note that the control field reaches 
significantly deeper than any part of the central field (despite its shorter exposure time) 
because of its much lower background light. In practice, we approximate each of the curves 
in Fig. ^ by the two-parameter Pritchet interpolation function ([Fleming et al. 1995|) 



where m (the "limiting magnitude" in a formal sense) is the magnitude at which / = 0.5, 
and a measures the steepness of decline of f(m) near m . The fitted values of (a, m ) for 
each annular zone and the control field are listed in columns (3) and (4) of Table 3. Column 
(5) gives the total area of each zone lying within the boundaries of the central field. 



1 



1 



a(m — itlq) 



(1) 
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The artificial-star tests were also used to evaluate the internal photometric precision, 
with the results shown in Fig. Our limiting magnitudes m as listed in Table 3 correspond 
roughly to the levels at which the measurement uncertainty reaches ~ 0.25 mag. 



3. THE LUMINOSITY DISTRIBUTION: RADIAL TRENDS? 

McLaughlin (1995) has used our HRCam data to analyze the radial structure of the 
inner GCS and to discuss its implications. Here, we investigate the characteristics of the 
GCLF (the luminosity function of the globular clusters) in the GCS core and compare it 
with previously published data for the outer halo. 

Early studies of the GCLF in the Milky Way and M31, along with the first photometry 
of the brighter clusters in the Virgo elliptical galaxies, gave rise to the now-standard notion 
that the GCLF is roughly Gaussian in number of clusters per unit magnitude (e.g., |Hanes| 
1977| ; |Harris fc Racine 1979| ; [Harris 1991| ). On purely empirical grounds, the near-uniformity 



of the GCLF peak or "turnover" absolute magnitude has encouraged its use as a standard 
candle QJacoby et al. 1992| [Harris 1996b] ; |Whitmore Hgg [Kavelaars fc Hanes 1997| ). The 



first photometric study that penetrated unequivocally past the turnover point in any 
galaxies beyond the Local Group was by Harris et al. (1991) for the Virgo ellipticals NGC 
4472 and 4649. By now there are published studies for several other large galaxies that 
reach clearly fainter than the GCS turnover level and directly verify the near-Gaussian 
shape of the GCLF (Fleming et al. 1995| ; |Whitmore et al. 1995| ; Forbes 1996a,b).Q 



In each of the four zones of our central field, the GCLF is defined by the number 
of objects remaining after (a) removal of nonstellar images by image moment analysis, 
(b) correction for detection incompleteness /, and (c) subtraction of the background LF, 
also corrected for incompleteness and normalized to the area of the zone. The results are 
summarized in Table 4, in 0.4-mag bins for each ring. The first column gives the central R 
magnitude of the bin, and the next four columns list the residual number of objects in each 
zone, fully corrected for incompleteness and background subtraction. For comparison, the 
last column gives the background LF (also corrected for incompleteness). From the relative 
zone areas listed in the last row, it can be seen that the background corrections are much 



3 It should be realized that the Gaussian distribution is only a convenient fitting function and has no 
astrophysical justification. Other equally simple analytic forms have been found that match at least as well 



(e.g. the ts function of Becker 1992). As Harris (1991) points out, if the GCLF were truly Gaussian in 
form, then so too would be the luminosity-weighted luminosity function (LWLF). However, the LWLFs are 
decidedly non- Gaussian ( McLaughlin 1994 : McLaughlin fc Pudritz 1996| ) . 
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less than 10% of the GCS population in every zone and magnitude range except for the two 
faintest bins in Ring 4. Note also that we stopped the calculations in each zone when the 
completeness / dropped below ~ 0.4, thus the effective limiting magnitude increases with 
radius. 

A potential concern is that our derived LFs might be contaminated at the faint end by 
false detections from pure noise, particularly in Ring 1 where the background light is high. 
To measure this effect explicitly, we inverted the image of the central field (i.e., we reversed 
the sign of all the pixel values on the sky-subtracted image) and then carried out exactly 
the same photometric analysis on this inverted picture as on the original one (see Harris 
et al. 1991 for a more extensive description of this procedure). The detected objects on this 
inverted frame are all pure noise by definition. The vast majority are either fainter than 
our adopted limiting magnitude in that radial zone, or are rejected as nonstellar in shape, 
and the remaining corrections to the LFs are quite small, amounting on average to about 4 
objects in the faintest bin in each zone. These corrections have been included in the final 
totals given in Table 4. 

The GCLFs for the separate rings are plotted in Fig. |5|. For rings 2, 3, and 4 the GCLF 
turnover (peak point) expected at R ~ 23.2 has clearly been reached and passed. The 
question of immediate interest is whether or not any systematic differences in the features 
of the GCLF show up with galactocentric distance. Standard Kolmogorov-Smirnov (K-S) 
two-sample tests show that rings 2, 3, and 4 are not significantly different from each other. 



This result extends that of [McLaughlin et al. (1994)1 , who found no variation in the GCLF 



over the radial range 1'.2 < r < 6'8. The same conclusion can also be drawn from attempts 
to fit Gaussian curves to the GCLFs, in which we find that the turnover point Rq and 
dispersion o~(R) of the fitted Gaussians are virtually identical for these three outer zones. 
Combining rings 2 — 4 for R < 23.8, we find by weighted least squares that the best-fit 
Gaussian has Rq = 23.18 ± 0.25 and o = 1.39 ± 0.15. This fit is displayed in Fig. |^. Adding 
the mean color (V — R) = 0.50 (see below), we obtain Vo (turnover) ~ 23.7, in complete 
agreement with the values Vo = 23.72 ± 0.06, a = 1.40 ± 0.06 quoted by Whitmore et al. 
(1995) from an HST-based sample that reaches about half a magnitude deeper than ours. 

For the innermost zone (ring 1), our data reach their reliable limit just at R ~ 23.4, 
and thus we cannot make any stringent test of the GCLF turnover level except to conclude 
that R is no brighter there than in the outer zones. Although Whitmore et al. (1995) do 
not comment directly on this point, they also apparently did not find any strong changes 
in the turnover out to their radial limit of r(max) ~ 114". If we assume Rq — 23.2 for 
ring 1 as well, then the data give a weak hint that the dispersion may be narrower there, 
with a best-fit o ~ 0.8 ± 0.3. This result is indicated schematically in the lower panel of 
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Fig. |6]. However, the relatively small sample in ring 1 prevents us from strongly ruling out 
a higher a. A K-S test comparing ring 1 with either ring 4 or the sum of (2+3+4), over the 
magnitude range R < 23.4, indicates that they are different at the ~ 90% confidence level - 
again, suggestive but not definitive. 

An alternate way of displaying the data to compare more readily with theory is as 
dN/dL, the number of clusters per unit luminosity; we call this function the LDF (HP94)fj 
Theoretical formation models, as well as a large body of observational material for currently 
active star- and cluster-forming regions (e.g., HP94; [McLaughlin fc Pudritz 1996| ; |Elmegreen| 



fc Falgarone 1996] ; [Elmegreen fc Efremov 1997] ), suggest that the initial mass distribution 



function of the clusters should be closely approximated by a power-law form dN/ dM ~ M ~ 7 
where the mass spectral index 7 is typically in the range ~ 1.5 — 2.0. The LDF will have 
the same form as long as the mass-to-light ratio for individual clusters is independent of M. 

We display the LDF in Fig. [7|, for three distinct regions covering virtually the entire 
range of the M87 halo. Our innermost zone (ring 1) is shown as the upper set of points, 
covering the projected radial range ~ 0.7 — 1.4 kpc, while rings 2-4 (~ 1.4 — 10.3 kpc) 
are plotted as the middle set of points. Finally, for the LDF of the outer halo, we use the 
V— band GCLF data of [McLaughlin et al. (1994)| covering the radial range 2'.88 — 6f82 



(~ 12.5 — 29.8 kpc) and shown as the lower set of points. 

The principal dynamical mechanism of tidal shocks should remove clusters more 
efficiently at lower mass and thus progressively flatten the LDF over time - and, of course, 
more rapidly at smaller R gc . How large is this effect expected to be? The Murali & 
Weinberg (1997) simulations starting from a very plausible initial mass function indicate 
that 7 should decrease by ~ 0.3 over 10 Gyr of evolution in an M87-like core. In our data, 
we find that over the luminosity interval 4.7 < log(L/L ) < 6.2 (i.e., clusters more luminous 
than the GCLF turnover point), direct weighted fits give 7 = AlogiV/AlogL = (1.59 + 0.16) 
for the outer-halo zone, (1.65 ± 0.17) for the mid-halo, and (1.84 ± 0.48) for the inner zone. 
No significant differences emerge in any of the regions, and if anything, the innermost zone 
exhibits a steeper mass function. This latter result is not an artifact of increased background 
noise at faint magnitudes, since we explicitly measured the (small) corrections due to noise 
and subtracted them out (see above). f\ Our data, with their admitted limitations, therefore 



4 As noted by McLaughlin (1994) and McLaughlin & Pudritz (1996), the GCLF in its conventional form 
as number of clusters per unit magnitude is just the luminosity-weighted first moment of the LDF. 

5 The larger uncertainty on 7 for the innermost zone is obviously a result of its much smaller sample size. 
Unfortunately, the sample cannot be increased, since our data already include every cluster brighter than 



the turnover point and within r < 30" of the nucleus. The outer-halo sample of McLaughlin et al. (1994) also 
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show no clear evidence for the trends indicated by the models. 

Alternately, we might ask more specifically about the effect of dynamical evolution 
on the classic GCLF turnover point (plotted in its normal observational form as number 
per unit magnitude). The progressive removal of low- mass clusters in the inner regions, 
as noted above, should cause the turnover point to become brighter with time and the 
GCLF dispersion to become narrower. The most recent numerical simulations which deal 
with the GCLF (Capriotti & Hawley 1996, Vesperini 1997, Ostriker & Gnedin 1997) 
find that the GCLF shape tends to be roughly preserved as evolution continues, once it 
approaches its present Gaussian-like form, and that the evolution in the turnover point can 
be surprisingly modest. In particular, Ostriker & Gnedin (1997) predict that the cluster 
populations in the inner ~ 5 kpc of large galaxies like M31, M87, and the Milky Way 
will have turnover points that are brighter by typically (~ 0.3 ± 0.1) magnitude than the 
outer-halo populations in the same galaxies, along with narrower GCLF dispersions by 

typically 0.1 0.2 magnitude. These radial differences - both rather modest - are at 

least weakly consistent with our observations, which show similar turnover levels to within 
±0.2 in all of our zones, but a narrower dispersion in the innermost zone. We remark that 
Kavelaars & Hanes (1997) find a rather similar observational result for the Milky Way and 
M31: the inner-halo GCLFs are more sharply peaked, yet have turnover luminosities that 
are statistically indistinguishable from those in the outer halo. An earlier version of this 
same conclusion was found by [Armandroff (1989~)| , and the steady improvements in the 



Milky Way database have reinforced it. 

In summary, (i) the modern theoretical simulations predict rather modest effects on the 
GCLF with radius (at least for the brighter part of the GCLF that is most easily observed), 
and (ii) our data indicate trends with radius that are, if anything, even smaller than the 
model predictions. The destructive mechanisms described above are important primarily on 
clusters fainter than the nominal turnover (log L/L Q < 4.7, or 1O 5 M ). Thus to find 
more unequivocal evidence for dynamical effects, we will need to explore the GCLF to much 
fainter magnitude levels and at widely separated places throughout the M87 halo. Note, 
however, that the deepest available study ([Whitmore et al. 1995|) shows no clear evidence 



for radial variations in the GCLF for My < — 5.4, corresponding to M^3 x 10 M . 

A more extreme - and more contrived - possibility is that we are witnessing the 
effects of convergent evolution, in the sense that the initial cluster mass distribution was 



includes the great majority of the clusters brighter than the turnover. The only experimental prospect for 
drawing firmer conclusions on the radial dependence of 7 is to obtain similar observations for many more 
galaxies. 
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steeper in the core, and has just now reached a flatter slope which coincides with the more 
slowly evolving outer-halo LDF. Such an explanation seems unlikely to us particularly 
because it would have to apply equally well for many different types of galaxies (dwarfs, 
spirals, and ellipticals) within which the detailed destruction rates would differ. Also, the 
formation models of McLaughlin & Pudritz (1996) and Elmegreen & Efremov (1997), which 
successfully account for the power-law shape of the LDF, already give strong grounds for 
expecting that the initial LDF should be independent of galactocentric radius. 

Another notable feature of Fig. |7| is the very top end of the LDF (log L/L Q > 6, 
equivalent to M^2 x 10 6 M©), which cuts off more sharply for the middle and inner zones 
than for the outer halo; that is, the very most massive globular clusters seem progressively 
to disappear as we move inward. Is this feature showing us evidence of dynamical friction, 
which is expected to become important for exactly this mass range ( [lremaine et al. 1975| ; 
Capuzzo-Dolcctta fc Tesseri 1997] )? The recent simulations of Capriotti & Hawley (1996) for 
clusters in a large-galaxy potential predict that most of these massive clusters (10 6 — 10 7 M Q ) 
should survive if R gc > 20 kpc; about half will be destroyed near R gc ~ 8 kpc; and almost all 
will disappear within ~ 2 kpc. For M87, we find from Fig. |^ that there are 19 ± 8 clusters 
brighter than V ~ 20.25 (log L/L Q ~ 6) in the outer zone; 4 ± 3 clusters in the middle 
zone; and none in the inner zone. However, these numbers must be normalized to the same 
parent population: the total number of clusters over the range 20.25 < V < 23.25 is ~ 650 
in the outer zone, 350 in the middle zone, and 40 in the inner zone. Thus in proportion to 
these totals, we would have expected to see 10 ± 4 of the super-bright clusters in the middle 
zone if none were destroyed, but about half this many if the models are correct; we observe 
~ 4 in the real sample. In the inner zone, we would expect just one cluster if no dynamical 
destruction has occurred; we see none. Although these comparisons are only approximate 
ones (the radial binnings in the models are rough ones, the statistics are weak, and the 
observed numbers are for projected radial bins rather than true three-dimensional radii), 
the model predictions are consistent with the data. 

Another way to state this conclusion is that, for the radial regime (r> 1 kpc) that we 
can most easily observe, dynamical friction has had little effect on the cluster LDF. We 
therefore concur with Lauer & Kormendy (1986) and McLaughlin (1995) that the GCS 
core radius, which is ^ 4 kpc, is too large to have been created by the action of dynamical 
friction on a hypothetical GCS that initially followed the more centrally concentrated stellar 
light profile (unless the inner clusters all formed in severely radial, plunging orbits that took 
them well within the stellar bulge; see |(Japuzzo-Dolcetta fc Tesseri 1997], |Capuzzo-Dolcetta 



fc Vignola 1997|. However, we regard this as a highly speculative alternative lacking any 



direct evidence; see the papers cited for additional discussion). 
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Interestingly, to detect the effects of dynamical friction on the lower-mass parts of the 
LDF more directly (i.e., M< 1O 6 M ), the models cited above indicate that we would have 
to study the cluster mass distribution within a radial regime that penetrates into the core 
radius of the M87 stellar bulge light, i.e. to within r <0.5 kpc (or about 7"). Unfortunately, 
there appear to be no immediate prospects for doing this: at such small radii, it is currently 
impossible to find and measure the fainter clusters, and there are virtually no clusters 
present in the first place within this small projected area. In other galaxies, the GCS 
populations are much smaller than in M87, and the situation is therefore even worse on 
observational grounds. 

In conclusion, we favor the view that dynamical evolution has not had very dramatic 
effects on the shape of the GCLF for the mass range M> 1O 5 M , except possibly at the 
very top end (M ^ 2 x 10 6 M Q ) where dynamical friction may have cut off the most massive 
clusters within a few kpc of the nucleus. The best present observations are beginning to 
converge with current theoretical simulations. Still deeper measurement of the cluster 
population at several different places in the M87 halo would be extremely valuable, with 
the potential to place still more stringent limits on the dynamical models. 

4. COLOR AND METALLICITY DISTRIBUTIONS 

The shorter exposures of the central field in V, I (Table 1) were used to measure color 
indices for the brighter M87 clusters and thus to investigate their metallicity distribution. 
The color-magnitude diagrams from our data are shown in Fig. |[ From the luminosity 
function of the control field (Table 4), we estimate that just 10 objects brighter than 
R ~ 23 should be due to field contamination, so the points plotted in these diagrams are 
almost entirely M87 globular clusters. The overall mean color indices of the sample (for 155 
clusters brighter than R = 22) are (V — R) = 0.501 ± 0.008 with rms scatter o = 0.10 mag; 
(R- I) = 0.617 ± 0.011 with a = 0.13 mag; and (V - I) = 1.118 ± 0.008 with a = 0.17 
mag. The scatter due to photometric measurement uncertainties alone averages ±0.07 mag 
in each filter over this magnitude range. The intrinsic dispersion in cluster color should 
then be ~ 0.15 mag in (V — I), the index most sensitive of the three to metallicity. 

Recently, Elson & Santiago (1996a,b) and Whitmore et al. (1995) have obtained HST 
(V, I) photometry for samples of M87 globular clusters from three different regions in the 
halo. They find a distinctly bimodal color distribution - previously suspected to exist by 
Lee fc Geisler (1993)| from {C — T%) photometry with somewhat lower internal precision - 
with histogram peaks at (V — I) ~ 0.92 and 1.23. These colors correspond (see below) to 
[Fe/H] ~ —1.7 and 0.0. In addition, at high metallicities, all of the conventional photometric 



-14- 



indices may give systematically incorrect metallicity values (additional comments will 
be made on this point below). Regardless of the exact metallicities of each of the two 
groups, bimodal distributions of this type have the traditional interpretation that the 
halo formed in two major enrichment stages, involving comparable amounts of gas. Very 
similar bimodal distributions are clearly present in the GCSs of many giant ellipticals (e.g., 
Forbes et al. 1997| |Geisler et al. 1996| ; |Zepf et al. 19951) as well as in the Milky Way (Zinn 



1985; Armandroff 1989), but may be absent, or considerably less obvious, in several other 
ellipticals ([Ajhar et al. 1994] ; [Forbes et al. 199"6|) . Bimodality thus seems to be a common, 
but perhaps not universal, phenomenon (but see |C6te et al. 1998|) . 



Elson & Santiago (1996a,b) also claim the existence of a mild trend in mean metallicity 
with magnitude in the sense that the bluer (more metal-poor) clusters tend to be brighter. 
Similar trends are not obviously present in the much larger samples of clusters measured by 
Whitmore et al. (1995) or [Lee fc Geisler (1993) , and Forbes et al. (1997) do not find any 



such effect in another giant elliptical, NGC 5846. Careful inspection of Elson & Santiago's 
data (cf. Fig. 2 from their second paper) shows that their claimed trend in fact arises from 
an excess population of clusters in the "blue" group (V — / ~ 0.92) concentrated near 
V ~ 22.0. This extra little clump of objects is not near the top end of the GCLF, despite 
its appearance in the Elson-Santiago diagrams; it is actually almost two magnitudes fainter 
than the GCLF tip. 

Our sample of clusters is drawn from a parent population more than three times 
larger than the Elson-Santiago sample, so we should be able to make a stronger test for 
the reality of this feature. In Fig. ||, the mean (V — I) color is plotted as a function of 
magnitude in 0.5— mag bins, extending to the faint limit of our data. Our conclusions from 
this exercise are that (a) the very brightest clusters (based on just 7 objects brighter than 
R = 20) tend to fall in the "red" (more metal-rich) group; and (b) for R > 20, there is 
no statistically significant change in mean color with magnitude. Nevertheless, we detect 
the Elson-Santiago feature as a slightly bluer than average color (by 0.05 mag) at just 
the magnitude level that they noticed (R ~ 21.5, V ~ 22.0). In conclusion, we verify the 
existence of this mildly anomalous feature, but we do not find that it affects the overall 
structure of the GCLF in any important way. Nevertheless, their study raises the possibility 
of interesting fine structure in the color-magnitude diagram which should be pursued with 
higher precision photometry of a larger sample. 

Finally, we use our (V — I) measurements to investigate the overall metallicity gradient 
in the M87 halo. Our material, covering the inner halo, can be combined with that of |Lee| 
fc Geisler (1993)| for the outer halo to give us a comprehensive view of the entire galaxy. To 



convert (V — I) into metallicity, we assume E(V — I) = 0.03 flBurstein fc Heiles 1984|) for 
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M87 and use the conversion relation 

(V - J) = 0. 18 [Fe/H] + 1 .20 (2) 

which we derive from the recent database of colors and metallicities for the Milky Way 
clusters ([Harris 1996a|) ; it is nearly identical with earlier calibrations (Couture et al. 1990; 
Kissler-Patig et al. 1997|) , and is based essentially on the [Zinn ( 1985 )| metallicity scale for 



the Milky Way clusters. We note in passing that these photometrically derived metallicities 
depend critically on the assumption that this simple linear relation is valid even at 
metallicities [Fe/H] ~ 0.0 and higher, even though there are no Milky Way globulars with 
accurately measured metallicities this high. This simple assumption may not be correct and 
that the extrapolation of this relation may seriously overestimate the true metallicities for 
[Fe/H] > - 0.5 (see |Harris et al7~m% jCarretta fc Gratton 1997|) . 



Our results for the mean metallicities are listed in Table 5, for seven radial bins with 
~ 20 clusters per bin: column (1) gives the mean radius of the clusters in the bin, column 

(2) the number of objects falling within the broad range 0.8 < (V — I) < 1.5, columns 

(3) and (4) the mean color and the rms scatter of the sample, and columns (5) and (6) 
the mean metallicity and dispersion. These points are plotted in Fig. |TD|, in which the 
error bars shown are the standard deviations of the mean metallicities directly from Table 
5. The zero point uncertainty of the [Fe/H] scale (i.e., the external error) is expected to 
be < ± 0.25 dex, from the ±0.03— mag uncertainty in the photometric zeropoint of our 
(V — /) scale (§2.1 above) and the similar uncertainty in the coefficients of the conversion 
relation Eq. (2). In Fig. [H], we also plot the [Lee fc Geisler ( 1993 )| metallicity measurements 
derived from their {C — T\) index. Encouragingly, in the radial range of overlap between 
our data and theirs (r = 60" — 100"), the mean metallicity values agree to well within ±0.2 
dex in [Fe/H], suggesting that the calibrations of both photometric metallicity scales are 
self-consistent to the level we expected. 

The total range in radius covered by Fig. |H] now extends from 10" to 500" (or 0.7 kpc 
to 36 kpc). It is apparent that the mean GCS metallicity decreases steadily outward from 
r ~ 60" to the outermost limits of the data, even though the range of cluster metallicities 
is comparably broad at any given radius (Lee & Geisler 1993; |Whitmore et al. 1995| ; Elson 



& Santiago 1996a,b). Although the scatter of the mean points does not permit any very 
precise determination of the slope, the rough relation 

[Fe/H] ~ -0.9 logr" + 1.15 (r > 60") (3) 

adequately describes its systematic change. An interpretation suggested by Geisler et al. 
(1996) from their extensive analysis of the NGC 4472 GCS, which is quite similar to that 
of M87, is that the bimodal metallicity distribution is present at all radii and that there 
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is little, if any, metallicity gradient within either of the two subsystems. However, the 
metal-richer ones form a much more centrally concentrated subsystem (as they do in the 
Milky Way), so their relative numbers increase strongly inward and thus create the observed 
shift in mean metallicity with radius. Notably, however, we find that within the GCS core 
radius of r c = 1' 4.3 kpc, the mean metallicity stays constant with radius. 



5. THE FORMATION OF SUPERGIANT ELLIPTICALS 

5.1. M87 and the Specific Frequency Problem 

The preceding discussion indicates that for M87, the inner-halo part of the GCS is 
simply an inward extension of the rest of the halo: the GCS in the core has a luminosity 
distribution function scarcely different from that of the outer halo, and the core metallicity 
distribution has a mean [Fe/H] ~ —0.3 which would result from its being drawn primarily 
from the metal-richer part of the bimodal distribution, along with a smaller admixture from 
the metal-poor component. All the available data are quite consistent with the picture that 
the internal characteristics of the M87 GCS exhibit no abrupt changes with radius. 

To this point, we have not yet brought in another important characteristic of the M87 
system, which is the specific frequency Sn ( Parr is fc van den Bergh 198 1| ; [Harris 199 1| ), or 



number of clusters per unit galaxy luminosity. The specific frequency is 

S N = N GC ■ 10°- 4 ( M £+ 15 ) (4) 

where My is the integrated magnitude of the host galaxy and Nqc is the total number 
of globular clusters in the galaxy. It has long been realized that Sn correlates strongly 
with environment: the highest- Sjv systems are found only in some cD-type giants at the 
centers of rich galaxy clusters ( [Harris et al. 19951 ; [Blakeslee 19~97| [hereafter B97]), while the 



sparsest globular cluster populations usually reside in "field" or small-group ellipticals in 
lower-density environments ( Harris 199 1| : |Wcst 1993 ). Thus to develop a larger picture of 



the origin of the GCS, we must also connect M87 to similar central-giant ellipticals in other 
clusters. We now turn to this somewhat broader issue. 

We first provide an improved calibration of Sn for M87, by combining our new data 
with the wider-field counts of McLaughlin et al. (1993, 1994) and Harris (1986). McLaughlin 
et al. find 3100 ± 130 clusters brighter than V = 24 over projected radii 1'97 < r < 9'.09, 
while the GCS radial profile of McLaughlin (1995) from our HRCam data gives 1300 ± 80 
clusters brighter than V = 24 and within r = lf97. For the GCLF parameters adopted 
above (cr = 1.4, V = 23.7), the globulars brighter than V = 24 account for 58% ± 6% of the 
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whole population; thus, the total number of M87 clusters with r < 9'09 = 39.7(D/15) kpc 
is 7590 ± 1050. Then from the M87 surface photometry of |de Vaucouleurs fc Nieto (1978)| , 



we find that the integrated magnitude of the galaxy interior to r = 9'.09 is My = —22.3 
[for (B — V) — 1.0 and (m — M)y = 31.0]. The metric specific frequency of M87, which is 
defined by B97 as the ratio of cluster numbers to galaxy light inside r = 40 kpc, is then 

5^°(M87) ~ (7590 ± 1050) x i °- 4 (-22.3+i5) = 9.1 ± 1.3 . (5) 

The number density of globular clusters in M87 falls off more gradually with r 
than does the intensity of background galaxian light ( [Harris fc Smith 1976| ; [Harris 1986 



McLaughlin et al. 1993|) , so the specific frequency must increase with galactocentric radius 
and the global Sn of the entire galaxy necessarily exceeds Sf}. To estimate the former 
quantity, we adopt r = 25' = 109(D/15) kpc as the "edge" of the galaxy (which has a total 
magnitude of My = —22.43; de Vaucouleurs & Nieto 1978) and of the GCS. The observed 
surface density profile of the outer cluster system can be taken from the wide-field counts of 
Harris (1986); these show that the cluster population beyond 9' makes up (1.048 ± 0.081) 
times the sum from r = 2' to 9'. The total number of clusters over all radii and magnitudes 
is then 13200 ± 1500, and hence 

S N (M87) = 14.1 ±1.6 . (6) 
The main source of uncertainty is the behavior of a c \ in the outermost halo of M87. 

The global specific frequencies of other elliptical galaxies in Virgo are typically of order 
Sn ~ 5 ( [Harris 1991| ). A valuable template for a 'normal' galaxy in this respect is M49 
= NGC 4472, which is comparable to M87 in brightness and size, but which has a global 
Sty = 6.0 ±1.5 and a metric ~ 4 ± 1 (from the data of [Harris 19SB]). Over all types 



of environments, Sn for E galaxies ranges from a minimum below Sn ~ 1 in some field 
ellipticals up to ^ 15 for the richest cD's. The underlying cause for this amazingly wide 
range among otherwise similar galaxies has been rather a mystery since the phenomenon 
first became clear two decades ago (e.g., [van den Bergh 1977| ; [Harris fc Petrie 1978| ; |Harri£ 



van den Bergh 1981|). 



5.2. Scaling Relations for Brightest Cluster Galaxies 

The search for an empirical set of necessary and sufficient conditions for high Sn has 
been a long one. While the highest— Sn galaxies often have a cD morphology, and are 
always the central members of their parent clusters or subclusters, there are many other cD 
and brightest cluster galaxies (BCGs) with quite normal GCS populations (e.g., McLaughlin 
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et al. 1994; Harris et al. 1995). However, the recent work of B97 and Blakcslee et al. (1997) 



(hereafter BTM97) considerably enlarges the database of GCS properties for these galaxies. 
Blakeslee et al. show that Sff in the BCGs is correlated with several indicators of the total 
cluster mass, such as the velocity dispersion of the galaxies and the X-ray luminosity of the 
hot gas residing in the cluster potential well. Moreover, it is now apparent from the BTM97 
study that BCGs cover a smooth continuum of specific frequency values, from normal levels 
up to the most extreme known. 

First we collect the available data for these BCGs. BTM97 give the "metric" Sff 
values for BCGs in 19 Abell clusters, measured by their surface brightness fluctuation 
technique. In addition, Harris et al. (1995) give global Sn values for other BCGs measured 
by direct resolution of the GCS population. From seven galaxies measured through both 
techniques (NGC 1399, 3842, 4486, 4874, 4889, 6166, and 7768), we derive a mean ratio 
(Sn/Sx) = 1.3 ± 0.2 which can be used to convert approximately between these two 
quantities in cases where only one has been measured. In Table 6, we summarize the 
updated parameters for 11 centrally dominant giant E galaxies in clusters not listed by 
BTM97. Here, columns (1) and (2) give the galaxy name and host cluster; column (3) the 
total luminosity of the BCG (assuming H Q = 80 km s _1 Mpc -1 ); columns (4) and (5) the 
global and metric specific frequencies; column (6) the velocity dispersion a c i of the galaxies 
in the cluster; and column (7) the X-ray temperature of the hot ICM gas in the cluster. 
Reference sources are listed in parentheses in columns (6-7). Adding them to the Blakeslee 
sample, we now have a total of 30 BCGs spanning the entire range from quite sparse groups 
(Leo, Centaurus) all the way up to the richest types of Abell clusters.^ 

The basic difference between the BCGs and other E galaxies is shown graphically in 
Fig. 11. Here, the total population of the GCS (calculated from the global Sn and the 
galaxy luminosity) is plotted against My. A total of 82 elliptical galaxies are now included 
in this plot; data for the non-BCG systems are taken from the compilations of [Harris & 



[Harris (1998)| and Durrell et al. (1996) with a few more recent updates from the literature. 
Despite the significant scatter, it is clear that to first order, N GC is directly proportional 
to galaxy luminosity (that is, Sn ~ constant) over a range of more than ~ 10 4 in L ga i. 
However, for the BCGs alone, Sn is itself a function of galaxy size, increasing steadily 
toward larger and more massive systems. 

As is clearly shown by B97 and BTM97, BCG specific frequency also increases 



6 In a few cases, the BCG is actually not the optically brightest galaxy in the cluster (Virgo and Coma 
are well known examples in which another cluster elliptical is slightly more luminous than the central one). 
Here, we are using the term "BCG" specifically to mean the centrally located giant E galaxy in the cluster. 
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with total cluster mass (represented either by the velocity dispersion a c i or the ICM gas 
temperature T x ). In Fig. [12| we show Sn versus T x for our augmented sample of BCGs, 
and Fig. |T3] is a similar plot of Sn and Nqc versus a c i. The overall dispersion of points in 
both graphs shown here - as well as the several others shown in BTM97 - is larger than for 
the Blakeslee sample by itself. Nevertheless, the general trend remains that more populous 
GCSs and higher specific frequencies are found in galaxy clusters with systematically higher 
mass, deeper potential wells, and hotter and more massive ICMs. Probably the clearest 
representation of these trends appears in the plots against a d . As did BTM97, we adopt 
a c i as the primary indicator of the size of the cluster, and using this, we can derive several 
simple empirical scaling relations for specific frequency, total GCS population, and BCG 
luminosity. Least-squares solutions against cluster velocity dispersion give the following: 

\ og S N = (0.71 ± 0.26) \oga d - (1.11 ± 0.68) , (7) 

log Nqc = (1.58 ± 0.32) \oga d - (0.42 ± 0.86) , (8) 

log (L gal /L Q ) = (0.86 ± 0.19) logo* + (8.55 ± 0.51) , (9) 

where a c \ is in km s _1 . In other words, we have approximate scaling laws Nqc ~ o" 1 ' 6 ) 
L ~ cr ' 9 , and Sn ~ N/L ~ a 0,7 . (We ascribe no particular significance to the fact that 
these correlations all pass through (0,0); at the low- mass end, all these relations require 
downward extrapolations which are not constrained by our data.) These relations stand 
independently of any of the theoretical hypotheses that will be discussed below. 

Much of the scatter in all these correlations may simply be due to the sizable 
measurement uncertainties in the individual Sn values. Only two objects stand out as 
genuinely anomalous: these are NGC 1399 (a high— Sn but low-luminosity cD in the 
relatively low-mass Fornax cluster), and NGC 1275 (the extremely luminous but Iow—Sn 
central cD in the massive Perseus cluster). The populous GCS in NGC 1399 has been 
measured several times with highly consistent results ( Hanes fc Harris 1986 : [Bridges et al 



1991| ; [Kissler-Patig et al. 1997| , among others), and the galaxy's luminosity profile is well 



enough established (e.g. [Killeen fc Bicknell 1988| ) to permit accurate calibration of either Sff 



or Sn (global). NGC 1399 is a challenging anomaly in almost any formation scenario. In the 
case of NGC 1275, the central galaxy is clearly extremely luminous even after accounting 
for the recent star formation present there, and the GCS is just as clearly unexceptional 
( |Kaisler et al. 1996| ). NGC 1275 might plausibly be explained by supposing that the Perseus 



cluster is a rare instance where the central cD built up from an initially rather normal E 
galaxy through amalgamation with many neighboring galaxies, as discussed by Kaisler et al. 
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5.3. Evolutionary Scenarios 

How did these intriguing BCG systems arise? The ideas that have developed in the 
literature can be grouped into three main approaches: (1) the high-Sjv galaxies acquired 
large numbers of "extra" globular clusters in various possible ways long after their initial 
formation stage; or (2) the high SV is a form of optical illusion, in which a population 
of intergalactic globular clusters that simply occupy the potential well of the cluster as 
a whole lies at the same location as the BCG; or (3) the GCS populations were largely 
built in at birth, with later minor adjustments. In the next three sections, we discuss each 
of these options in turn. A key to keep in mind is that any proposed solution to the 
problem must explain not only the sheer total number of clusters present in these galaxies; 
it must also address several other measurable GCS characteristics including their spatial 
distribution, velocity distribution function, and metallicity distribution function. As yet, 
no complete model has been developed which can quantitatively deal with all these aspects 
of the problem. To help evaluate the various options, we first comment on the evolutionary 
scenarios. 

Among the proposed mechanisms for supplementing the GCS populations of dominant 
cluster galaxies relatively late in their lives are the following: 

(1) Cooling flows: M87 and other cD galaxies have extended coronae holding large amounts 
of hot low-density gas, and it has been hypothesized that globular clusters somehow 
actively form out of the ongoing cooling flow from this gas, thus adding to the original 
cluster population ( |Fabian et al. 1984|) . However, this model fails on several grounds to 



predict the quantitative features of the observed GCSs in these galaxies, such as the cluster 
metallicity distribution; the complete lack of correlation between GCS population and 
cooling-flow amplitude; the near-complete absence of any young, massive star clusters in 
most cooling-flow galaxies; and the large velocity dispersion in the GCS which is unlike that 
expected from objects condensing out of a cooling flow (see |Grillmair et al. 1994b ; Harris 



et al. 1995; |Kaisler et al. 1996t [Bridges et al. 1996| ; |Holtzman et al. 1996Q 



(2) Stripping from other galaxies: the possibility that central BCGs might have acquired 
halo globulars from other, fully formed galaxies during close encounters was explored in 
a series of models by Muzzio and collaborators (see [Muzzio 1987| for a review). These 



investigations showed that cluster swapping may be responsible for minor rearrangements 
of clusters between galaxies, but is unable to produce the major increases in Sn (factors 
of 2 to 3) that we need for the central cD galaxies. A related possibility, in which tidal 
debris from systems throughout a galaxy cluster sinks onto a BCG at the bottom of the 
potential well (again, after the main epoch of galaxy formation), is discussed and rejected 
by |McLaughlin et al. (1994] . 
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Another version of this approach, pointed specifically at explaining the bimodal 
metallicity distributions seen in giant elliptical GCSs, has recently been developed by Cote 
etal (1998). They assume that a typical giant elliptical such as NGC 4472 started as a 'seed 
galaxy' of smaller size, and has since absorbed large numbers of smaller (mostly dwarf) 
galaxies. The metal-richer globular clusters in the final gE are thus mostly ones that formed 
with it originally, while the metal-poorer clusters are the ones acquired from the accreted 
dwarfs. By quantitative numerical simulations, Cote et al. show that a plausible range 
of GCS metallicity distributions can be accounted for this way. To account for a much 
higher-SAT BCG like M87, they assume further that the giant central galaxy accretes only 
the high-S^r envelopes of the smaller galaxies orbiting in the potential well of the galaxy 
cluster, without absorbing the central (cluster-poor) cores of those galaxies. The advantage 
of this general scheme is that it employs a process (small-galaxy accretion and harassment) 
which must surely happen at some level in almost any rich environment like Virgo, Fornax, 
or Coma. However, some related worries are that (i) to explain the entire GCS metallicity 
distribution function this way, it is necessary to assume that very large numbers of small 
galaxies are accreted (many hundreds for a case like NGC 4472, and even more for M87); 
(ii) if (typically) half the light of the gE is from smaller, accreted, and more metal-poor 
systems, they should not have the high (~ solar or higher) mean metallicities that they are 
observed to have; and (iii) an environment like Virgo should contain hundreds (perhaps 
thousands) of stripped cores of dwarf galaxies whose envelopes were accreted by the central 
BCG. These expectations run counter to the available observations. Nevertheless, their 
discussion correctly calls for additional investigation through detailed N-body simulations. 

(3) Mergers: E galaxies can be built out of the collisions and mergers of disk galaxies, 
and massive globular-like star clusters can form in the shocked gas during the merger 
(Schweizer 1987, 1997; [Kumai et al. 1993a|) . The newly formed clusters may end up in a 



more centrally concentrated and more metal-rich subsystem than the original (metal-poor, 
spatially extended) GCSs of the colliding galaxies ( |Ashman fc Zepf 1992| ; |Zepf fc Ashman 



|1993| ). Several galaxies in which this process is actively happening have now been identified 
(e.g., |Holtzman et al. 1992| ; |Whitmore et al. 1993b| ; |Whitmore fe Schweizer 1995| ; |Hortzmanj 



\et al. 1996 ; Miller et al. 1997 ). It is important to keep in mind that the specific frequency 



of the merged remnant may end up either larger or smaller than in the pre-merger galaxies, 
depending on factors such as the total amount of incoming gas, the collision geometry, and 
(most importantly) the efficiency of cluster formation relative to field-star formation. 

The merger mechanism is arguably quite a plausible way to build ellipticals that 
end up with low-to-normal Sn (Harris 1981, 1995; [Miller et al. 199?]) , since typical disk 



galaxies have relatively few halo clusters of their own, as well as modest supplies of gas. For 
example, isolated E Galaxies such as NGC 3557 ( Morbey fc McClure 19851) and NGC 720 
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( |Kisslcr-Patig et al. 1996 ), which have elongated shapes and very low specific frequencies 



in the range Sn<2, are prime candidates to be the results of disk-galaxy mergers that 
happened long ago. The more recent merger remnants such as NGC 3921 ( [Schweizer et al. 



1996 ) and NGC 7252 ( Miller et al. 1997 ) appear to be evolving into the same type of 



elliptical, with modest populations of globular clusters and Sjv ~ 2 — 3, as predicted by 
Harris (1995). However, applying the same model to giant ellipticals at the highest end of 
the Sn scale leads to much more serious difficulties. For such galaxies, thousands of massive 
star clusters would have to be formed during the mergers, requiring an enormous input 
supply of gas - so large, in fact, that the merging "galaxies" must in essence be pregalactic 
clouds ( |Harris 1995] ). In addition, some high-SAT galaxies do not have clearly bimodal GCS 
metallicity distributions; the highest-SAT systems do not have the lowest radial metallicity 
gradients; and the lower-metallicity clusters in these gE's are often not as metal-poor as 
those in spirals. All of these points - and others raised in the references cited above - speak 
against the merger model as a viable origin for the extreme BCGs ( |Gcisler et al. 199B] ; 
Forbes et a/.^997|) . 



In summary, these secondary-evolution processes appear able to produce measurable 
changes in Sat for ellipticals with cluster populations in the low or normal range. However, 
for M87 and the other supergiant cD galaxies with the highest specific frequencies, no viable 
production mechanisms have convincingly emerged from these approaches. 



5.4. Intergalactic Globular Clusters 



White (1987^ suggested that clusters of galaxies might host populations of globular 
clusters that are not currently bound to any one galaxy, but instead move freely in the 
potential well of the cluster as a whole. |West et al. (1995~)| have revived this idea, and 
cite what is essentially a version of Blakeslee's (1997) correlation between BCG specific 
frequency and galaxy cluster X-ray luminosity as possible evidence for the existence of 
these intergalactic globular clusters (IGCs). In their scenario, the Sat of a central BCG 
could be artificially enhanced by the superposition of IGCs — which would naturally have 
their largest space density in the very core of the cluster — on top of a normal (Sat ~ 5) 
population of globulars that actually are bound to the BCG. Indeed, given the existence 
of diffuse light in galaxy clusters (e.g. Coma: |Thuan fc Kormendy 1977|) , and the recent 
identification of candidate intergalactic stars in Fornax ( |Theuns fc Warren 1997] ) and Virgo 
( |Arnaboldi et al. 1996| ; |Ferguson et al. 1996| ; |Mendez et al. 19"9"7| ; |Ciardullo et al. 199^ ). 
it would be surprising if IGCs do not exist. However, the question whether a population 
of IGCs could single-handedly give rise to the high-SAr phenomenon is rather a different 
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matter. 



In our view, and as discussed by White (1987)1 , any IGCs should be found only in 



association with diffuse cluster light from intergalactic stars. That is, we consider it highly 
unlikely that globular clusters preferentially condensed out of the hot, diffuse gas that 
pervades galaxy clusters, leaving no other form of stellar material. Furthermore, to generate 
a high Sj\f in the central galaxy, this IGC component must itself have (as we will see 
quantitatively below) an extremely high and nonuniform Sn- A plausible origin for this 
component of an intracluster medium could be from the stripping of material from young, 
gas-rich galaxies during the collapse of a cluster (see also [Merritt 1984| ), at an early enough 



stage that many galaxies had formed globular clusters but had yet to produce most of their 
halo field stars. Most of the residual gas in these young galaxies would then have been 
heated to the virial temperature of the surrounding cluster and joined the ICM which, in 
large clusters, outweighs the stellar material in the galaxies ( Pavid et al. 1990 ; Arnaud 



et al. 1992] ; [Whitmore et al. 1993a|) . This option is worth investigating in some detail, since 



it stands potentially to address, at once, the issues of BCG specific frequency, diffuse light 
in galaxy clusters, and cD envelope formation. 

The first problem is to account for a correlation between cluster X-ray luminosity and 
the specific frequency of the intergalactic stellar material (and hence for the correlations of 
West et al. and B97). Let us define a globular cluster formation efficiency 



N GC 



M st ars + M gas 



(10) 



where M stars is the mass in the stellar component of the galaxy and M gas is the mass of the 
leftover gas unused for star formation. We can then write the presently observed specific 
frequency as 

N G c .._ (.. , M gas 



SN(X w^ = e 1 • • ( u > 

Wl stars \ m stars/ 

Suppose now that e is basically similar from place to place (see HP94; we will also use this 
assumption in the next section). Then, any galaxies in which star formation proceeded 
nearly to completion (M gas /M stars <C 1) would end up with essentially equal (and normal) 
Sjy set by the value of e. However, many slow-forming, late-type galaxies could well have 
been destroyed at a point in their lives when M gas was still large. The stars and globular 
clusters from these shredded gas-rich objects would also have spilled into the cluster at 
large, producing an Sn larger than normal by a factor ~ (1 + M gas /M stars ). In absolute 
terms, we could suppose that the collapse of more massive galaxy clusters engendered 
the destruction of larger numbers of gaseous systems, thus leaving behind larger ratios of 
gas mass to total stellar mass, as is observed (|David et al. 199"0| ; |Arnaud et al. 19*9^ ) . As 
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well, this effect could contribute to the growth of cD envelopes and the establishment of a 
correlation between cD envelope luminosity and the X-ray luminosity Lx of the surrounding 
cluster gas flSchombert 1988| ; [Lopez-Cruz et al. 1997[ ). 

Pursuing this option, let us now assume that a population of IGCs with specific 
frequency Sff , and an envelope composed of diffuse intracluster light L env , are 
superimposed on a BCG body with normal S l ^ dy . The apparent specific frequency of the 
combination will be 

qbody | / j it \ qIGC 
ntot _ J N ~r \Mervol ^bodyl^N /-, \ 

N ~ TZ7 77 ■ ^ ll > 

1 T J^env/ ^b 



Since Sffi c > S b ^ dy by hypothesis, the addition of the IGC component will increase S£°* 



L env /Lb dy increases weakly with the parent cluster Lx flSchombcrt 1988| ). Thus, Sffi would 



increase with Lx even if Sjfc were constant from cluster to cluster. However, Sffi in cD 
galaxies should then also correlate with the relative size of their cD envelope, L env / L^y. 
This is not an observed trend (for example, M87 itself has a rather modest envelope but one 
of the highest and best-measured specific frequencies; and other BCGs with very massive 
envelopes such as NGC 6166 have unexceptional Sn values; see the next section). 

The alternative is to let Sff also increase with cluster Lx (or mass, or velocity 
dispersion, or temperature). By Eq. (|TT|), this will occur if the formation of more massive 
galaxy clusters somehow involves the destruction of protogalaxies with systematically larger 
gas fractions M gas /M stars . In this case, the initial metallicity of the intracluster gas should 
be antz-correlated with Lx (see also Pavid et al. 1990] ). Even though some subsequent 



contamination due to galactic winds from the surviving protogalaxies must certainly have 
occurred (e.g., |Matteucci fc Gibson 1995|), such an anti-correlation does appear to exist 



( [Fabian et al. 1994]) . Thus, this scenario would allow IGCs to resolve at least some aspects 
of the high-S^ problem if we allow Sj^ c to increase with cluster mass. Unfortunately, this 
descriptive scenario cannot account quantitatively for other features of the GCS, as we will 
see next. 

Let us try to construct a model for M87 specifically, by building it from an M49-like 



body with a normal GCS, surrounded by a much more extended IGC envelope. In Fig. 14 



we model the observed 5-band surface photometry of M87 ( Carter fc Dixon 1978| ) as the 



sum of two components. One is the underlying body of the galaxy, with a conventional r 1//4 
law profile, 

^body _ i^ qi^ _|_ 2.658 (r/arcsec) 1//4 mag arcsec -2 . (13) 

The other is a faint envelope which, for the sake of our argument, is assumed to be part of 
the intracluster medium around M87: 

fji% v = 26.675 + 1.875 log [1 + (r/1200") 2 ] mag arcsec~ 2 . (14) 
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The envelope necessarily has a very large core radius of 24'.7 ~ 108(D/15) kpc. Beyond 
this core radius, the envelope light behaves roughly 5 . We note that Carter & 

Dixon (1978) find an especially bright and extensive envelope component by comparison 
with de Vaucouleurs Sz Nieto (1978) or other authors; thus, our two-component model is a 
"maximal-envelope" solution. 

Integration of Eqs. (13) and (14) shows that in the radial regime V<r<6'.5, over 
which Cohen & Ryzhov (1997) have recently obtained high-quality radial velocities for 205 
globular cluster candidates in M87, the envelope contributes to the total light of the galaxy 
in the ratio L env /Lb d y ~ 0.10 — 0.15. The total specific frequency of M87, averaged over 
this same interval, is roughly Sjv(M87) ~ 11, while that of M49 (which we shall take to 
represent the underlying S b ^ dy of M87) is ~ 4.7. Thus the fraction of the globulars around 
M87 that are actually bound to it would be 



Nbody _ S N f Lent- \ _ n ( -]i ,. 



l + ^0.4, l'<r<6f5. (15) 



tot D N \ lj body i 

The globulars bound to the galaxy will have a Gaussian velocity distribution with an 
intrinsic dispersion o v ~ 270 km s _1 characteristic of the stellar core of the galaxy (e.g. 
Sembach fc Tonry 1996|) , while the IGCs that are spatially superimposed will have a v ~ 570 



km s _1 typical of the early-type galaxies in Virgo (Binggeli et al. 1987] ). The combination 
of the two, in the proportions given by Eq. (fL5f) , can be reasonably well described by a 
Gaussian with o~ v ~ 400 km s _1 , which is not inconsistent with the data of Cohen & Ryzhov 
(1997). However, the most noticeable effect of any IGCs would show up in the wings of the 
velocity distribution: the Cohen-Ryzhov sample of 205 globulars built in these proportions 
should have ~ 5 clusters with v > 1000 km s -1 , whereas the actual sample has none. 

A still more stringent test comes from the decomposition of the M87 surface brightness, 
together with the observed specific frequency, into putative "body" and "envelope" (IGC) 
components. Suppose that the body of M87 has a specific frequency profile S b ^ dy identical 
to that of our 'template' galaxy M49 (shown in Fig. [l^, along with the total Sff observed 
for M87). By hypothesis, the difference between Sffi and S h ^ dy is due to the intergalactic 
envelope, which allows us to find Sffi c as a function of r. The results of this exercise are 
shown in Fig. |I~5|b as the solid dots. The open squares in the same figure show the Sff 



that is inferred if the "body" of M87 is assumed instead to have Sn = 5 throughout, with 
no dependence on galactocentric radius. 



The implications of this exercise are severe. In particular, we find that Sffi c oc 



-1±0.25 



between r ~ V.2 and r ~ 9' (~ 5 to 40 kpc). Because the cD envelope, or diffuse cluster 
light, has an essentially constant surface brightness at these radii (see Fig. [14] ), this means 



that Sj^ c is directly proportional to the projected number density o~igc of the supposed 
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intergalactic globulars; thus, we conclude that oigc ^ Rgc ±0 ' 25 - That is, we are forced to 
claim that the intergalactic clusters are significantly concentrated towards the center of 
the galaxy itself, with a core radius no larger than ~ V = 4.3 kpc (that is, the same core 
radius as the M87 GCS itself). Such behavior is unacceptable for any collection of truly 
intergalactic objects, and completely contradictory to our original hypothesis. It is also 
worth noting that the deduced absolute value of Sffi reaches a level in the core (> 100) 
that is an order of magnitude larger than in any known galaxy. 

We might try to get around the strong radial dependence of the IGC cluster component 
by assuming an Sff c that is independent of galactocentric radius, and then use this 
and the observed Sffi and S b ^ dy to solve for the ratio L env /L body . The observed total 
surface-brightness profile of M87 can then be used to find the envelope intensity I env as a 
function of r. However, this rather desperate argument results in a diffuse light component 
that varies essentially as I env oc r~ l down to ~ V distances from the center of M87 - which 
again implies that this "intergalactic" material must in fact be bound to the galaxy. Finally, 
we note also that if we had used the de Vaucouleurs & Nieto (1978) photometry instead, 
the solution for the envelope (Fig. FJ) would have been considerably fainter, S\ff c even 



larger, and thus even less favorable to the scenario. 

It is important to emphasize more descriptively why the IGC model fails. The key is 
that the "specific frequency problem" embodied by M87 is a local problem everywhere in 
the halo as well as a global one. M87 has an anomalously large cluster population at all 
galactocentric radii r and not just in its cD envelope. Thus, for example, we cannot turn 
M49 into M87 simply by adding large numbers of extra clusters to its outskirts; they must 
be added in the same proportions all the way in to the inner halo and core of the galaxy 
as well ( |Harris 1986| ). This requirement then forces the postulated IGC profile to have a 



central concentration resembling the underlying galaxy, as well as an inordinately high 
specific frequency. There appears to be no natural way to do this through the IGC model, 
because a realistic IGC component, or a population of clusters tidally stripped and accreted 
from other galaxies, would preferentially boost the cluster population in the outer halo of 
the galaxy while leaving the inner halo relatively unchanged. 

The argument we have just outlined may, however, not apply to all high-Sjv BCGs. 
Some other BCGs have relatively much larger cD envelopes and thus, probably, stronger 
contributions from the central cluster potential well. The intriguing case of NGC 1399, the 
BCG in the Fornax cluster, may represent one such situation. Recent kinematic studies of 
its GCS QKissler-Patig 1998j |Minniti et al. 1998| ) show that for radii > 20 kpc the velocity 



dispersion of its globular clusters matches the surrounding Fornax galaxies rather than the 
NGC 1399 halo light (a difference of a factor of two in o~(v r )). We are clearly just beginning 
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to understand how these rather unusual galaxies were built. 



5.5. A Formation Scenario 

The preceding arguments favor the view that the essential features of the M87 GCS, 
and the other high-Sjv cD systems like it, were built in at very early times. Somewhat 
by default, this "initial conditions" route has usually relied on the supposition that the 
efficiency of globular cluster formation e as defined in Eq. ( |10D could have differed by about 
one order of magnitude ab initio among E, cD, dwarf, and spiral protogalaxies, thus directly 
creating the large range in Sn that we now see (e.g., Harris 1981, 1991; [van den Bergh 1984| ; 



HP94). It is possible that the amount of gas cloud shocking and turbulence could have 
directly governed e, driving it up to higher values in denser or more massive protogalactic 
environments where collision speeds were higher. Unfortunately, to date this has been a 
somewhat ad hoc view with little observational evidence that could be brought to bear one 
way or the other, and so it, too, has been a less than satisfactory option. 

The recent work of B97 and BTM97 opens the door to an alternate approach to this 
problem. Their data (see above) imply more clearly than before (see West 1993] and [Kumai 



\et al. 1993] for similar earlier suggestions) that Sn increases with the overall mass density 



of the protogalactic environment. B97 proposes essentially that the high-S 1 ^ galaxies do 
not have excess numbers of clusters; instead, they have artificially low total luminosities 
relative to their globular cluster populations. In this picture, the "missing light" is in the 
hot intracluster gas that never formed stars. 

In such a scenario, then it is natural to ask why this gas was unused for star formation, 
and how it ended up where it did. The earliest ~ 1 Gyr history of an E galaxy is expected 
to be heavily influenced by the SNII supernova rate from the first generation of massive 
stars, and the subsequent development of a galactic wind ( [Vlathews fc Baker 1971| ; [Larson] 



1974| ; see |Loewenstein fc Mushotzky 1996| or Gibson 1994, 1996, 1997 for recent overviews). 



Under the right conditions a large fraction of the original protogalactic gas can be ejected 
outward by the wind to join the ICM. It seems highly likely to us that galactic winds would 
have influenced the GCS specific frequency, and by different degrees in different galaxies. 
Our reason for proposing this is based primarily on the timing of the GCS formation epoch: 
a large protogalaxy should accumulate from many smaller (10 8 — 10 9 M Q ) dwarf-sized gas 
clouds (supergiant molecular clouds or SGMCs), which provide the right host environments 
for globular cluster production ( |5earle fc Zinn 1978j ; Larson 1990b, 1993; HP94). In this 



picture, we postulate that each SGMC contains several protocluster cores, which are the 
much denser sites where star formation can take place at the necessary > 50% efficiency 



-28- 



to produce a bound star cluster. Since the protocluster cores are expected to build up 
over timescales ~ 10 8 y out of the gas in SGMCs (HP94; [McLaughlin fe Pudritz 199(j| ; 



Elmegreen fe Efremov 1997]), it is likely that they would form very early in the evolution 



of the protogalaxy, before the lower-density majority of the gas began star formation in 
earnest. Thus, the globular clusters should already be present before the wind reaches its 
peak at a few x 10 8 y and interrupts subsequent field-star formation. 

In essence, we suggest that for normal E galaxies with specific frequencies Sn ~Z 5, 
the early galactic wind did not eject large amounts of the original gas supply. Today, 
these galaxies have GCS specific frequencies fairly close to the original cluster formation 
efficiency. However, in systems like M87, we propose that the galactic wind was stronger, 
ejecting much of the original gas mass outward into the dark-matter potential well of the 
surrounding cluster, and leaving behind a GCS with an artificially high Sn- 

Is this idea a realistic one for gas within such massive potential wells as we find in the 
BCGs? (Normally, we would expect the less massive protogalaxies to be able to eject higher 
proportions of their initial gas.) |Gibson (1997)| provides a comprehensive discussion of the 



way that various parameters affect the SNII production and mass ejection rate. Effective 
ways to increase the amount of mass loss are (a) to increase the star formation efficiency 
factor v = M+/M tot , which would boost the early supernova rate; or (b) to decrease the 
rate at which the expanding supernova remnants can cool radiatively and lose energy; or 
(c) to adopt an IMF flatter than the normal Salpeter (x = 1.35) slope, thus increasing the 
number of Type II SNe per unit mass. Other studies have suggested that an IMF slope 
x ~ 1 is also appropriate for matching the observed Fe and O abundances in the ICM gas 
(e.g. [Loewenstein fe Mushotzky 1996| ; |Arimoto fe Yoshii 1987| ; [Renzini et al. 1993| ; Pavid 



\et al. 199T| , among others). We are, however, looking for a way to generate a difference in 
the mass loss rate (and hence Sn) in a BCG like M87, compared with normal systems like 
M49. Given that these galaxies have rather similar compositions and structures, we see no 
compelling reasons to suggest arbitrary differences in the SNR cooling rates, or the IMF. 
Thus of the various options mentioned above, the most plausible appears to us to be the 
star formation rate v. Our scenario therefore boils down to the suggestion that BCGs had 
an extraordinarily high v at early times. 

Another important point to consider is that the simple one- zone wind models mentioned 
above are an oversimplification for large protogalaxies, which were likely to have been highly 
clumpy. Thus the parametrizations of mass loss rates for monolithic models such as the 
ones reviewed by |Gibson (1997)| should not be taken too literally. For a central supergiant 



like M87, which is at the gravitational center of its surroundings, the protogalaxy would 
have been the site of a considerable amount of gas infall and chaotic, violent motion which 
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could have encouraged early, rapid star formation and a more massive SNII-driven wind. 
This picture is consistent with the considerable evidence that many large E galaxies almost 
certainly formed at epochs predating z ~ 3, and some probably began major star formation 



at z>5 (e.g., Larson 199Ua; Maoz 199U; 


Turner 1991; jWhitmore et al. 1993a 


; Loewenstein & 


Mushotzky 1996; 


Mushotzky & Loewenstein 1997; 


Stcidcl et al. 1996; |Giavalisco et al. 1996; 



Pender et al. 1996] ; [Ellis et al. 19971 , to name a few; see also HP94 for a timing argument in 
support of a redshift z> 5 for the first epoch of globular cluster formation). 

We emphasize again that we are not suggesting this approach as the way for a typical 
galaxy (elliptical or spiral) to form. The BCGs are, clearly, rather unusual entities, and the 
formation processes during the much higher-density regime of z ^ 5 are entirely likely to 
have been different from those at z < 2 which gave rise to the main population of galaxies 
that we now see. A fully realistic model of a BCG-type protogalaxy during the epoch of 
globular cluster formation and galactic wind outflow, in all its undoubted complexity, is far 
beyond the scope of this paper (or indeed any current models). However, we can roughly 
quantify the effect on the specific frequency as follows: from Eqs. (f|), (|10D, and ([11]), we 
obtain 



S 



N 



85.5 x 10 ! 



M. 



gas 



(16) 



L J V V M s tars 

where [MjV)y is the mass-to-light ratio for the stellar material in an E galaxy; we adopt 
(M/L) v = 8 flFaber fc Gallagher 19791 ; IBinney fc Tremaine 19871) ■ M gas is the amount of 
gas assumed to have been driven out by the galactic wind and unused for star formation. 
With the explicit assumption e ~ const, Eq. ([II]) can be rewritten as 
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where S N represents the fiducial specific frequency for a galaxy in which little or no gas is 
ejected and the star formation runs to completion (M gas -C M stars ). Lacking a full physical 
theory, we must estimate S N on observational grounds: 

(a) The mean specific frequency for 14 large ellipticals in Virgo, Fornax, and the NGC 5846 
group (from data in Harris & Harris 1997; see also [Harris 1991[ [Kissler-Patig et al. 1997| ) is 
(Sn) = 5.4 ± 0.3, excluding the cD's M87 and NGC 1399. Since even the normal ellipticals 
have probably lost some of their original gas to early winds, this value should give an upper 
limit to S%. An extreme lower limit can be set at S% ~ 2, which is the mean for the lowest 
specific frequencies found in a range of E and dE galaxies in sparse groups. A straight mean 
over all the non-BCG galaxies, weighted inversely as the observational uncertainty in Sn, 
gives (S N ) = 3.3 ± 0.2. 

(b) The mean specific frequency for 13 dwarf ellipticals more luminous than My ~ — 15 
( Purrell et al. 1996|) is (Sn) = 4.2 ± 0.5. These same dE's are small enough to have lost 
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typically > 20% of their initial gas in a wind, according to published models (e.g. Gibson 
1997) J] However, they may also have later destroyed a significant fraction of their original 
number of globular clusters by dynamical friction, since the GCS in a dwarf elliptical is 
quite centrally concentrated (see the discussion of Durrell et al). These two effects have 
opposite and comparable influences on Sn, and it is not immediately clear which will be 
the more important over a Hubble time. 

Taking the preceding arguments to bracket our estimate of S N , we adopt S N = 3.5 ± 1. 
(Similar 'baseline' values for S N were reached by West et al. 1995 and BTM97, though for 
different reasons.) From Eq. fll6|), we find that the corresponding fiducial cluster formation 
efficiency is e ~ 5.8 x IO^Mq 1 , or one globular cluster per 1.7 x 10 8 M Q of gas. For a 
mean cluster mass ~ 3 x 1O 5 M (from HP94; [Harris 19964 |Meylan fc Heggie 1997| ), the 



corresponding specific mass ratio as defined by HP94 is Sm — Macs/M stars = 1.8 x 10~ 3 . 
This result agrees extremely well with the observations from contemporary star-forming 
regions in the Milky Way and M31, in which the mass of gas seen to be forming clusters in 
giant molecular clouds is about 2 x 10~ 3 of the total GMC mass (HP94; McLaughlin 1997| ). 

The value of S% indicates that only 0.2% of the stellar mass in a typical E galaxy 
today is in the form of bound star clusters. To obtain a case like M87 with Sn ~ 14, 
we must then require (M gas /M sta rs) — 3, i.e. that 75% of the original gas supply in the 
protogalaxy was ejected to rejoin the ICM. More generally, if M gas equals the amount 
expelled by the wind, then the fraction f M of the initial protogalactic gas that is lost 
is f M = M gas / {M gas + Mstars) = 1 - (S%/S N ). This quantity, for our BCG sample, is 
shown in Fig. pl| For comparison, the expected dependence of / on a derived from the 
scaling relations Eqs. (7-9) above is shown by the dashed line. A rough linear relation 
I'm ~ (cr c ;/1200kms _1 ) as shown in the graph matches equally well, given the observational 
scatter. 

From this graph, we see again that to save the assumption e ~ const, we have to 
claim that quite a substantial fraction of the initial gas mass was expelled from the 
biggest systems. In terms of the total gas mass in the ICM, this is not at all an excessive 
requirement, since large Abell-type clusters typically contain ~ 10 13 — 10 14 M Q of diffuse 
gas, most of which is likely to be primordial material never processed through galaxies (e.g. 
David et al. 1991] ; |Whitmore et al. 1993a| ). We recognize, however, that this requirement 



imposes very strong demands on the early SNII rate which will have to be put to the test 



7 Among dwarf E galaxies with globular cluster systems of their own, the lowest-luminosity dE's have the 



highest specific frequencies (Durrell et al. 1996), the most extreme cases being the Fornax and Sagittarius 
dwarfs with Sn ~ 25. These tiny systems are likely to be ones which expelled most of their original gas in 
a single early burst flDckcl fc Silk 1986| ). 



in more advanced modelling. 

Lastly, we explore the dependence of total GCS size on the total core mass in the 
surrounding cluster including the dark matter. For an isothermal potential well, the central 
density is given by p c = 9a 2 j 'AnGr 2 where r c is the projected core radius that we directly 
observe. The cluster core radii are known empirically (BTM97; |Schaeffer et al. 1993| ) to 
scale approximately as r c ~ 104 kpc (a/500km s^ 1 ) ' 6 . Thus we obtain 



M core is of course dominated by the dark matter, with smaller contributions from the 
visible galaxy and the X-ray gas. Since the core radii are typically of order ~ 100 kpc, 
p c defined in this way comfortably encloses the bulk of the central BCG and provides a 
reasonable representation of the overall mass density there. In Fig. [17] we show the ratio 
Mcore/Ncc, the total mass per unit globular cluster, plotted against velocity dispersion. The 
graph indicates that the globular cluster formation efficiency in an absolute sense (number 
per unit total mass, M core = M gas + M stars + M dark ) is lower in the more massive BCGs. 
This result is consistent with our suggestion that the galactic wind may have been earlier 
and more effective in the more massive environments, possibly interrupting even the GCS 
formation before it could run to completion. 

This conclusion about the absolute efficiency of globular cluster formation is not the 
same as in BTM97. This apparent disagreement is due primarily to different definitions of 
terms. They defined an efficiency ratio using a projected core mass M c over a cylindrical 
volume, derived from the surface density S c rather than p c . In addition, they calculated 
M c for a fixed radius of 40 kpc in all clusters. We believe it is more appropriate to use a 
total core mass calculated from a true three-dimensional mass density, in order to compare 
with the total globular cluster population in the entire BCG. For the sake of more direct 
comparison, we show in the lower panel of Fig. [T7| the ratio M 100 /Ngc, where M 100 is the 
total mass within a fixed (spherical, not cylindrical) volume of radius of 100 kpc, which 
roughly represents the total optical extent of these galaxies. In the sense defined by M 100 , 
the absolute efficiency increases with a c i, though the trend is less obvious than in the upper 
panel. It is not clear to us which of the two representations is preferable; in either graph, a 
typical mass ratio is about 3 x 1O 9 M per globular cluster. Comparing this value with the 
mean e _1 ~2 x 10 8 M Q (gas mass per globular cluster) derived above, we see that the dark 
matter makes up > 90% of the total mass even within the cluster core. 




(18) 



and the total mass contained within the core is 




(19) 
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The scenario we have outlined here explores the implications of a single important 
assumption: that the conversion rate e of protogalactic gas into globular clusters was 
statistically uniform in different environments. Pursuing this assumption leads to the 
hypothesis that BCGs ejected a much higher fraction of their initial gas mass than did more 
normal E galaxies. We recognize that this discussion is highly speculative, and is supported 
only in an indirect sense by the GCS scaling relations which are just now beginning to 
emerge from the data. 

The obvious alternative to this admittedly risky conclusion is simply to assume that e 
is strongly dependent on environment, varying by up to an order of magnitude in different 
E galaxies. As yet, there are no definitive arguments either to rule out such a view or to 
clearly support it. The galactic wind hypothesis is offered as an alternate approach, which 
does has the advantage of connecting to a known major event in the early history of a large 
E galaxy. 

6. SUMMARY 

We have conducted a photometric study of the inner-halo globular cluster system in 
the Virgo giant M87, and have used the characteristics of its GCS as the starting point for 
a discussion of the formation of central cluster galaxies like it. A summary of our principal 
findings is as follows: 

(1) The luminosity distribution of the globular clusters in M87 (in its form either as the 
GCLF or the LDF) is virtually independent of location in the halo for any projected radius 
r > 1 kpc. We suggest that the main effects of dynamical evolution (tidal shocking and 
evaporation) may be limited to clusters less massive than M C ;<1O 5 M , below the faint 
limits of the present data. Alternatively, cluster evolution has proceeded in such a way 
as to preserve the overall shape of the GCLF as a function of time, as the recent model 
simulations are beginning to show. Dynamical friction may have removed the very most 
massive clusters (M>2 x 10 6 M Q ) in the innermost few kiloparsecs. 

(2) Within ~ 4 kpc of the galaxy center, the mean metallicity of the GCS is uniform with 
radius at the metal-rich level of [Fe/H] ~ —0.3. At larger radii, the mean cluster metallicity 
declines as (Z/Z Q ) ~ r -°- 9 . We find that cluster metallicity is largely uncorrected with 
cluster luminosity. 

(3) M87 holds the paradigmatic high-specific-frequency globular cluster system. We 
examine carefully the hypothesis that its many thousands of "extra" globular clusters are 
due to an intergalactic (IGC) population in the Virgo potential well, and reject it. The 
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dominant quantitative problem with the IGC model is that it requires the intracluster 
globulars to have a highly unrealistic spatial distribution, inconsistent with the assumptions 
of the model. We conclude that the vast majority of the globular clusters within M87 
belong to it, and were formed with it in situ. 



(4) Adding the recent results of Blakeslee (1997) and |Blakeslee et al. (1997)| to other data 
from the literature, we reinforce the the result that GCSs in "brightest cluster galaxies" like 
M87 become continuously more populous (higher specific frequency) with parent galaxy 
size, and with the total mass of the surrounding galaxy cluster. However, the absolute 
efficiency of GCS formation (measured as the number of globular clusters per unit mass in 
the cluster core) decreases with increasing core mass. We use the available observations for 
30 BCGs to construct empirical scaling relations describing the globular cluster population 
and specific frequency as a function of cluster mass (measured by the velocity dispersion o c i 
of the galaxies). 

(5) We suggest that globular cluster specific frequency in BCGs may have been strongly 
influenced by the early galactic wind driven by the first round of Type II supernovae 
throughout the protogalaxy. In particular, we speculate that BCGs differed from normal 
ellipticals by generating a higher early star formation rate. The considerably stronger 
resulting wind interrupted subsequent star formation, drove out a large fraction of the initial 
gas mass, and left behind a galaxy with a higher than average ratio of globular clusters to 
field stars. If this hypothesis proves to be invalid, the most obvious alternative is to require 
that the GCS formation efficiency per unit gas mass (e) was higher in BCGs than in normal 
ellipticals, and in general becomes progressively higher in larger- density environments. 
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Figure Captions 

Fig. 1. — R— band image of the central field in M87. This frame is a median of 26 individual 
500-sec exposures in R taken with the CFHT and HRCam at prime focus (see Sec. 2.2 of 
the text). North is approximately at the upper right corner, East at lower right; see Fig. 2 
for a finder chart with an exact orientation map. The large black dot along the left edge of 
the field is the hole enclosing the bright guide star for the tip-tilt correction of the camera. 
The smooth isophotal contours of the galaxy have been subtracted through an ellipse-fitting 
model (see Sec. 3). 

Fig. 2. — Distribution of starlike objects in the M87 central field that are brighter than the 
GCLF turnover (i?< 23.2). Larger symbols correspond to brighter objects; compare with the 
direct image of Fig. 1. Fiducial directions are marked on the borders, with the galaxy center 
at (0,0). The concentric circles have successive radii of 9'/2, 19'/7, 39'/3, and 59'/0, defining the 
four radial zones used in the GCLF analysis. Note the small sector directly above the galaxy 
center with no measured objects: this sector contains the nuclear jet and was excluded from 
any analysis. 

Fig. 3. — Completeness of image detection in each annular zone, as a function of R 
magnitude. These curves are the empirical results of artificial-star tests as described in 
Sec. 2.2. The five zones are defined in Table 4: Ring 1 is the innermost zone and B denotes 
the remote background field. 

Fig. 4. — Internal precision of the R— band photometry, as derived from the artificial-star 
tests. Here <j{R) is rms scatter in the measured magnitudes of the recovered stars. The data 
are plotted (unsmoothed) in 0.5-mag steps, with the five zones labelled as in the previous 
figure. 

Fig. 5. — GCLFs for the four zones defined in Table 4. In each bin, the number of globular 
clusters per 0.4-mag bin, fully corrected for detection incompleteness and background 
subtraction, is plotted against R magnitude. 

Fig. 6. — (a) Upper panel: The GCLF for the three zones 2-4 combined. A Gaussian curve 
with peak at R = 23.2 and dispersion a = 1.4 mag as derived from Whitmore et al. (1995) 
and the present data is superimposed on the datapoints. (b) Lower panel: The GCLF for 
the innermost ring. Two Gaussian curves with the same peak but different dispersions are 
compared: the narrower curve with a = 0.75 provides the formal best fit, but the wider one 
which matches the outer zones cannot be strongly ruled out. 

Fig. 7. — Luminosity distribution function (LDF) for the M87 globular clusters in three 
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regions of the halo. Here dN/dL represents the number of clusters per unit luminosity, 
plotted against cluster luminosity in Solar units (L/L ); the three sets of points have been 
shifted vertically by arbitrary amounts for clarity. Solid circles are the data from ring 1 
(projected radii 0.7 — 1.4 kpc); open circles are from rings 2 — 4 (1.4 — 10.3 kpc); and solid 
triangles are outer-halo GCS data (r > 10 kpc) from |McLaughlin et al. (19941 Note the 



similarity of the LDF slope in all three regimes over most of the luminosity range, coupled 
with a progressive disappearance of the most luminous clusters (log (L/L Q ) > 6.2) going 
inward toward smaller radii. 

Fig. 8. — Color-magnitude distribution for the measured objects in the M87 central field; R 
is plotted against (V — R), (R — I), and (V — I). In the last panel, the two peaks in the 
bimodal color distribution noted by |Elson fc Santiago (1996b~)| and Whitmore et al. (1995) 
are indicated by dashed lines. Most of our objects belong to the redder, more metal-rich 
group since they are in the inner halo; see text. 

Fig. 9. — Mean (V — I) colors of the globular clusters as a function of magnitude. The 
dashed arrow at R ~ 21.5 marks the magnitude near which [Elson fc Santiago (1996b)] found 



an extra "clump" of clusters in the blue part of the bimodal color distribution. 

Fig. 10. — Metallicity gradient for the globular cluster system in M87, combining the inner- 
halo data from this study (solid symbols) and the outer-halo data from [Lee fc Geisler (1993" 



(open symbols). Note that within one core radius r c ~ 1', there is no change in the mean 
cluster metallicity; at larger radii, there is a steady decrease in mean metallicity out to the 
limits of the observations, scaling roughly as Z ~ r -°- 9 . 

Fig. 11. — Total number of globular clusters plotted against galaxy luminosity, for elliptical 
galaxies. Open circles indicate giant ellipticals which are the centrally dominant members 
of their clusters (BCGs). Solid dots indicate all other E-type galaxies, including dwarfs. 
The solid line is a direct proportionality Nqc ~ L ga i for a specific frequency S% = 3.5. See 
Sec. 5.4 of the text. 

Fig. 12. — Specific frequency Sn plotted against the X-ray temperature of the ICM gas. Open 



symbols are BCGs taken from Table 6. Solid symbols are BCGs taken from the |Blakcslcc 



\et al. (1997)| sample, where we have converted their metric Sff to the global specific frequency 



via Sn = 1.3Sn (see Sec. 5.4). 

Fig. 13. — (a) Upper panel: BCG specific frequency plotted against the velocity dispersion 
o~ci of the galaxies in the cluster. Solid symbols are the BTM97 sample, and open symbols 
are taken from Table 6 as above. The solid line shows the scaling relation from Eq. (6). (b) 
Lower panel: Total cluster population N GC ~ L ■ plotted against act. Symbols are as 
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above. The solid line shows the scaling relation from Eq. (17). 

Fig. 14. — Radial profile of the M87 halo light, deconvolved into 'body' and 'envelope' 
components as described in Sec. 5.4 of the text. Solid dots represent the B— magnitude 
surface intensity profile of Carter & Dixon (1978); open symbols represent the B 
measurements of |de Vaucouleurs fc Nieto (1978)[ The extended envelope component, defined 
to fit the Carter-Dixon data, is assumed to be intracluster light in the overall potential well 
of the Virgo cluster; see text. 

Fig. 15. — (a) Upper panel: Specific frequency profiles for the globular cluster systems of M87 
(solid symbols) and M49 (open symbols), as described in Sec. 5.4 of the text. A model for the 
M87 GCS is assumed here in which the main body of the galaxy has a normal M49-like specific 
frequency, while the intergalactic globular clusters from the extended envelope raise the total 
specific frequency up to the observed value Sffi. Both sets of data are taken from Fig. 21 of 
[McLaughlin et al. ( 1994 )| , and have been properly scaled to reflect a Virgo distance modulus 



of (m — M) v = 31.0 and the same GCLF turnover and dispersion as adopted in Sec. 3. (b) 
Lower panel: Deduced specific frequency profile for the IGC (envelope) component, obtained 
from the difference between the top and bottom curves in panel (a). The assumed IGCs are 
required to have an extremely high overall Sjv, and to be concentrated toward the galaxy 
nucleus. 

Fig. 16. — Mass fraction fu = 1 — (S%/Sn) plotted against o c i. According to the galactic 
wind hypothesis discussed in Sec. 5.5, Jm is the fraction of the initial protogalactic gas lost to 
the ICM, fu = M gas /(Mg as + M stars ). The dashed line indicates the predicted dependence 
of / on a from the scaling relations in Eqs. (7-9), while the solid line shows the more 
approximate linear relation defined in the text. 

Fig. 17. — (a) Upper panel: Ratio of total core mass in the cluster of galaxies to globular 
cluster population, plotted against velocity dispersion. See Sec. 5.5. (b) Lower panel: Ratio 
of "metric" core mass (within 100 kpc radius) to Nqc, plotted against velocity dispersion. 



